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QQ . Abstract 

We provide a justification with rigorous error estimates showing that the leading term in 
weakly nonlinear geometric optics expansions of highly oscillatory reflecting pulses is close to 
Ph ' the uniquely determined exact solution for small wavelengths e. Pulses reflecting off fixed non- 

■^^ I characteristic boundaries are considered under the assumption that the underlying boundary 

(-H ' problem is uniformly spectrally stable in the sense of Kreiss. There are two respects in which 

these results make rigorous the formal treatment of pulses in Majda and Artola |16j , and Hunter, 
Majda and Rosales [TO]. First, we give a rigorous construction of leading pulse profiles in prob- 
lems where pulses traveling with many distinct group velocities are, unavoidably, present; and 
second, we provide a rigorous error analysis which yields a rate of convergence of approximate to 
exact solutions as £ — > 0. Unlike wavetrains, interacting pulses do not produce resonances that 
^ ■ affect leading order profiles. However, our error analysis shows the importance of estimating 

pulse interactions in the construction and estimation of correctors. Our results apply to a gen- 
eral class of systems that includes quasilinear problems like the compressible Euler equations; 
iy^ I moreover, the same methods yield a stability result for uniformly stable Euler shocks perturbed 

l' ■ by highly oscillatory pulses. 
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1 Introduction 

We study highly oscillatory pulse solutions for a general class of hyperbolic equations that includes 
quasilinear systems like the compressible Euler equations. Our main objective is to construct leading 
order weakly nonlinear geometric optics expansions of the solutions (which are valuable because, for 
example, they exhibit important qualitative properties), and to rigorously justify such expansions, 
that is, to show that they are close in a precise sense to true exact solutions. 

A single pulse colliding with a fixed noncharacteristic boundary in an A^ x A^ hyperbolic system 
will generally give rise to a family of reflected pulses traveling with several distinct group veloci- 
ties. We study this situation when the underlying boundary problem is assumed to be uniformly 
spectrally stable in the sense of Kreiss. A formal treatment of this problem was given in Majda- 
Artola [16], building on an earlier treatment of nonlinear geometric optics for pulses in free space 
in Hunter-Majda-Rosales [10]. In the papers [161 HOj - systems of nonlinear equations for leading 
order profiles were derived, but their solvability was not discussed. Morever, the questions of the 
existence of exact solutions on a fixed time interval independent of the wavelength of oscillations 
(or pulse width) e, and of the relation between exact and approximate solutions, were not studied 
there. In this paper, we give a rigorous construction of leading pulse profiles in problems where 
pulses traveling with many distinct group velocities are, unavoidably, present. In addition, we 
construct exact solutions on a fixed time interval independent of e, and provide a rigorous error 
analysis which yields a rate of convergence of approximate to exact solutions as e — )■ 0. 

Rigorous treatments of the short-time propagation of a single pulse in free space were given in 
Alterman-Rauch [1| and Gues-Rauch [SjH . The methods (e.g., conormal estimates in [UlS], high- 
order approximate solutions in [8]) used in the constructions of exact solutions and in the error 
analyses of these papers do not readily extend to problems involving many pulses with distinct 
group velocities. The method we use here to construct exact solutions and justify leading term 
expansions involves replacing the original system (jl.ip with an associated singular system (11. 3p 
involving coefficients of order - and a new unknown Ue{x, 60)10 Exact solutions U^ to the singular 
system yield exact solutions to the original system by a substitution 



Ue{x) = Ue 

where <l)o{x') = x' • /3 is the "boundary phase" as in (jl.ip . Both the singular system and the system 
of profile equations satisfied by the leading profile U^{x, OQ,S,d) are solved by Picard iteration. 



^The paper 8 considered "fronts" as well as pulses. 

^The singular system approach was used in [2j in their study of a single pulse on diffractive time scales. 



The error analysis is based on "simultaneous Picard iteration", a method first used in the 
study of geometric optics for wavetrains in free space in [11]. The idea is to show that for every 
n, the n-th profile iterate Z^°'"(x,6'o, ^) converges as e — )• in an appropriate sense to the n- 
th exact iterate U^{x,6o), and to conclude therefrom that U^{x,6o, ^) is close to C/e(x,0o) for e 
small. Unlike wavetrains, interacting pulses do not produce resonances that affect leading order 
profiles. However, our error analysis shows the importance of estimating pulse interactions in the 
construction and estimation of correctors. Another key tool in the error analysis, discussed further 
in section 11.41 is the machinery of moment-zero approximations developed in section 14.11 Our use 
of these approximations was inspired by the "low- frequency cutoff' argument of • 

The main novelties of this paper are: 

1) We give a rigorous treatment of pulses refiecting off boundaries; earlier rigorous work con- 
cerned pulses in free space. 

2) We provide methods for handling many pulses traveling with distinct group velocities; in 
particular, we show that although pulse interactions do not produce new pulses at leading order, 
pulse interactions must be estimated in the construction of correctors and in the error analysis. We 
distinguish in the estimates between "transversal" and "nontransversal" pulse interactions. 

3) In contrast to the treatment of uniformly stable refiecting wavetrains in [6] , we are able here 
to give a rate of convergence of approximate to exact solutions as wavelength e — )• Oo 

1.1 Exact solutions and singular systems 

In order to study geometric optics for nonlinear problems with highly oscillatory solutions it 
is important first to settle the question of whether exact solutions exist on a fixed time interval 
independent of the wavelength (e in the notation below). A powerful method for studying this 
problem, introduced in [11] for initial value problems and extended to boundary problems in |24j . 
is to replace the original system with an associated singular system. 

On M^ = {x = {x',Xd) = {t,y,Xd) = {t,x") : Xd > 0}, consider the N x N quasilinear 
hyperbolic boundary problem: 

d 
^Aj{Ve)d^^Ve = f{Ve) 

(1-1) / , x' ■ 13 

b{ve)\xd=o = go + eG ix', 

Vg = uq in t < 0, 

where xq = t is time, G{x',9o) G C°°(M'^ x M^,MP) decays to zero as \6o\ — ^ oo, with supp G C 
{xo > 0}, and the boundary frequency /3 G M^ \ {ofl Here the coefficients Aj e C°°(M^,M^'), 
/ E C°°(M^,R^), and b G C°°(M^,MP). 

Looking for v^ as a perturbation v^ = uq + eue of a constant state uq such that /(uq) = 0, 



^In the case of wavetrains there was an "arithmetic obstacle" to obtaining a rate of convergence in the error analysis; 
namely, the generation of many noncharacteristic, but "almost characteristic", phases by nonlinear interactions. 
Because pulses interact weakly, that obstacle is absent in the problem studied here. 

^Wavetrains instead of pulses are obtained by taking G{x' , 6o) to be periodic in ^o- 



6(iio) = 50; we obtain for u^ the system (with slightly different Aj^s) 

d 

(a) P{eue,dx)us := '^Aj{eus) d^^Us = T{eue)ue on x^ > 

i=o 

(1-2) / x' ■ 13" 
(6) B{eUe)Ue\xa=o = G[x', 



(c) Ue = in t < 0, 

where i?(i^) is a C°° p x N real matrix defined by 

6(iio + eug) = b{uo) + B{eue)eue 

and J^ is defined similarly. We assume that the boundary {x^ = 0} is noncharacteristic, that 
is, Afi{Q) is invertible. The other key assumptions, explained in section [1.21 are that P{v,dx) 
is hyperbolic with characteristics of constant multiplicity for f in a neighborhood of the origin 
(Assumption 1 1 . 1| and that {P{0,dx),B{0)) is uniformly stable in the sense of Kreiss (Assumption 

[rnD. 

For any fixed Eq > 0, the standard theory of hyperbolic boundary problems (see e.g., [H [12]) 
yields solutions of (jl.2p on a fixed time interval [0, Teg] independent of e > Eq- However, since 
Sobolev norms of the boundary data blow up as e — )• 0, the standard theory yields solutions u^ 
of (jl.2p only on time intervals [0, T^] that shrink to zero as e — t- 0. In section [21 exact (and 
necessarily unique) solutions to ()1.2p of the form Ue^x) = Ue{x, 2—^) are constructed on a time 
interval independent of e G (0, Eq] for eo sufficiently small, where Us{x,Oq) satisfies the singular 
system derived by substituting U£{x, ^^) into (|1.2p : 

d rf-l 

(1.3) ^^^ ^^'^ 

B{eU,){U,)\x,=o = G{x',6o), 

U, = Omt<0. 



As explained in [23j, the study of singular systems is greatly complicated by the presence of 
a boundary. Even if one assumes that the matrices Aj are symmetric (as we do not here), there 
is no way to obtain an L^ estimate uniform in e by a simple integration by parts because of the 
boundary terms that arisqj- The blow-up examples of |^3j show that, at least in the wavetrain case, 
for certain boundary frequencies /3 it is impossible to estimate solutions of (jl.3p uniformly with 
respect to e in C{xd, H'^{x', 9o)) norms, or indeed in any norm that dominates the L°° nornu. We 
do not know if analogous blow-up examples exist in the pulse case, but it is clear that the proofs 
of this paper do not apply when f3 lies in the glancing set (Definition II. Sp . 

In [7] a class of singular pseudodifferential operators, acting on functions U{x',6q) decaying in 
^0 and having the form 

(1.4) Ps{Dx',eo)U := f e^'^'+^^^'^p feF(x', ^o),^' + — ,7) t/(C', A;) d^'dA;, 7 > 1, 

^The class of symmetric problems with maximal strictly dissipative boundary conditions provides an exception to 
this statement, but that class is too restrictive for some important applications; for example, the boundary problem 
that arises in the study of multi-D shocks does not lie in this class. 

®The problem occurs only for /3 in the glancing set fPefinition II. 3p . as the examples of [23j together with the 
results of [24] show. 



was introduced to deal with these difficulties. Observe that after multiplication by A^ i^Us) and 
setting Aj := A^^Aj, F := A'^^I', (L3\i becomes 



(1-5) = d.,Ue + A (eUe,d^,, + ^"j Ue = F{eUe)Us, 

B{£Ue){Us)\^,=0 = G{x',eo), 

Ue = Omt< 0, 

where A ( eU^, dx' H — ) is a (differential) operator that can be expressed in the form ()1.4p . Kreiss- 

type symmetrizers rs(-Dx',eo) i^ the singular calculus can be constructed for the system p.Sp as 
in [23] under the assumptions given below. With these one can prove Lp'{xci-,H^{x' ,6o)) estimates 
uniform in e for the linearization of (ll.5p . The main difference with |24j is that we use here a 
singular pseudodifferential calculus that is specially constructed for pulses in [7]. In the pulse case 
9q lies in an unbounded set and the exact profile [/^(xj^o) has continuous Fourier spectrum. The 
analysis of [23] relied on a singular calculus for wavetrains. In that case 6q lies in S^ and Ue{x^ Oq) 
has discrete Fourier spectrum, a fact that was used in several places for proving symbolic calculus 
rules in [23]- The results of the pulse calculus needed here are recalled in Appendix 1X1 

To progress beyond L'^{xd, H^{x' , Oq)) estimates and control L°° norms, the boundary frequency 
/3 must be restricted to lie in the complement of the glancing set (Definition [L3]). With this extra as- 
sumption we are able to use the pulse calculus to block-diagonalize the operator A ( sUe, dx' H — ^ J 
and thereby prove estimates uniform with respect to e in the spaces 

(1.6) E^ = C{xd, mix', eo)) n L'^ixd, H^+\x', Oq)). 

These spaces are algebras and are contained in L°° for s > -^. For large enough s, as determined 
by the requirements of the calculus, existence of solutions to ()A.6p in E^ on a time interval [0, T] 
independent of e G (0,eo] follows by Picard iteration (see Theorem II. 12p . 

1.2 Assumptions and main results 

Before continuing with an overview of the strategies for constructing profiles and for showing 
that approximate solutions are close to exact solutions, we now give a precise statement of our 
assumptions and main results. 

We make the following hyperbolicity assumption on the system (II. 2p : 



Assumption 1.1. The matrix Aq = I. For an open neighborhood U of E M , there exists an 
integer q > 1, some real functions Ai, . . . , Ag that are C°° on U x W^ \ {0} and homogeneous of 
degree 1 and analytic in ^, and there exist some positive integers vi, . . . ,Vq such that: 



det 



a q 



j=l ~ fc=l 



for u ^ U and ^ = {Cit ■ ■ i^d) G K'^ \ {0}. Moreover the eigenvalues Ai(n,^), . . . , Aq(n,^) are 
semi-simple (their algebraic multiplicity equals their geometric multiplicity) and satisfy \i{u,£,) < 
■■■ < Xq{u,0 for alluGU, ^eR'^X {0}. 



We restrict our analysis to noncharacteristic boundaries and therefore make the fohowing: 

Assumption 1.2. For u ^ U the matrix Ad{u) is invertible and the matrix B{u) has maxim,al 
rank, its rank p being equal to the number of positive eigenvalues of Ad{u) (counted with their 
multiplicity). 

In the normal modes analysis for the linearization of (jl.2p at G W, one first performs a Laplace 
transform in the time variable t and a Fourier transform in the tangential space variables y. We 
let T — i7 G C and ry G W^~^ denote the dual variables of t and y. We introduce the symbol 

A{O:=-iA-\Q){{T-i^)I + Y,V,A0)\ , C:=(r-i7,r?)GCxM^-i. 

For future use, we also define the following sets of frequencies: 

E := {(r -ij,n)eCx M-^-i \ (0,0) : 7 > o} , S := {c G H : r^ + 7^ + \rj\^ = l} , 

Ho := {(r, 7?) G M X R'^-^ \ (0, 0)} = H n {7 = 0} , So := S n Hq . 

Henceforth we suppress the u in Xk{u,^) when it is evaluated at n = and write Afc(0,.^) = Xk{0- 
Two key objects in our analysis are the hyperbolic region and the glancing set that are defined as 
follows: 

Definition 1.3. • The hyperbolic region % is the set of all {T,r]) G Hq such that the matrix 
A{t, rf) is diagonalizable with purely imaginary eigenvalues. 

• Let G denote the set of all (r, ,^) G M x M*^ such that ^ 7^ and there exists an integer 
k G {!,...,(/} satisfying: 

r + Xk{0 = ^iO=0. 

If vr(G) denotes the projection of G on the d first coordinates (in other words 7r(r, ^) = 
(r,^i, . . . ,£,d-i) for all (t,^)), the glancing set Q is Q := ir{G) C Hq. 

We recall the following result that is due to Kreiss |r2] in the strictly hyperbolic case (when all 
integers Vj in Assumption [TTT] equal 1) and to Metivier |19j in our more general framework: 

Proposition 1.4 ([12,,19j). Let Assumptions \L i\ and \l.S\ be satisfied. Then for all C € 'E\ Hq, the 
matrix A{C) has no purely imaginary eigenvalue and its stable subspace E,^{() has dimension pU| 
Furthermore, W defines an analytic vector bundle over H \ Sq that can be extended as a continuous 
vector bundle over S. 

For all {T,r]) G Hq, we let E*(t, r/) denote the continuous extension of E^ to the point (t, 77). The 
analysis in [191 shows that away from the glancing set Q C Hq, IE*((^) depends analytically on i^, 
and the hyperbolic region Ti does not contain any glancing point. 
Next we define the hyperbolic operator 

d 

i=i 
and recall the definition of uniform stability [121 '4]: 



^The stable subspace is the direct sum of the generahzed eigenspaces associated to eigenvalues with negative real 
part. 



Definition 1.5. The problem (jl.2p is said to be uniformly stable atu = if the linearized operators 
{L{dx),B(0)) at u = are such that 

B(0) : E*(r — i'j, rf) — )• C^ is an isomorphism for all (r — ij,rj) € S. 

Assumption 1.6. The problem (II. 2p is uniformly stable at u = 0. 



It is clear that uniform stability at u = implies uniform stability at nearby states (and therefore 
at all u gU up to restricting U). Thus, there is a slight redundancy in Assumptions 11.21 and 1 1 . 61 as 
far as the rank of -8(0) is concerned. 

Boundary and interior phases. We consider a planar real phase (pQ defined on the boundary: 

(1.7) Mt^y) ■=zt + v-y, {z,v)eEo. 

As follows from earlier works (e.g. |16j). oscillations on the boundary associated with the phase (pQ 
give rise to oscillations in the interior associated with some planar phases (pm- These phases are 
characteristic for the hyperbolic operator L{dx) and their trace on the boundary equals (J)q. For 
now we make the following: 

Assumption 1.7. The phase (pQ defined by (II. 7p satisfies {T,r]) G Ti. 

Thanks to Assumption 11.71 we know that the matrix A{t, rj) is diagonalizable with purely imaginary 
eigenvalues. These eigenvalues are denoted icoi, . . . ,iujMi where the w^'s are real and pairwise 
distinct. The Wm's are the roots (and all the roots are real) of the dispersion relation: 



d-l 



det 



rI + Y.r,A^{0)+u:Ad{0) 



-J 



0. 



To each root Um there corresponds a unique integer fcm S {1, • • • , g} such that r + Xk^{r],ujm) = 0. 
We can then define the following reajfl phases and their associated group velocities: 

(1.8) V?n = 1,... ,M, (pm{x) := (poit,y)+uJmXd, v^ := VAA..„(r/,a;m) . 

Let us observe that each group velocity v^ is either incoming or outgoing with respect to the space 
domain M!^: the last coordinate of v^ is nonzero. This property holds because (r, 77) does not 
belong to the glancing set G. We can therefore adopt the following classification: 

Definition 1.8. The phase (pm is said to be incoming if the group velocity v^ is incoming (that is, 
d^d^krr,i^,^m) > 0), and outgoing if the group velocity v^ is outgoing ('%Afc„(^,a;m) < 0). 

In all that follows, we let I denote the set of indices m E {1, . . . , Af } such that <pm is an incoming 
phase, and O denote the set of indices m S {1, . . . , M} such that (pm is an outgoing phase. If p > 1, 
then I is nonempty, while if p < A^ — 1, O is nonempty (this follows from Lemma [1.91 below). 

Main results. We will use the notation: 

d 

L(r,e):=rI + J^e,-4,(0), 

i=i 
/5=(r,r/), x' = it,y), M^') = P ' x' . 



*If (r, rj) does not belong to the hyperbolic region H, some of the phases 0™, may be complex, see e.g. |21II23|[T^ 
1181 19]. Moreover, glancing phases introduce a new scale ^/e as well as boundary layers. 



For each phase (j)m, dcpm denotes the differential of the function (pm with respect to its argument 
X = {t,y,Xd)- It follows from Assumption 11.11 that the eigenspace of A{(3) associated with the 
eigenvalue icOj^ coincides with the kernel of L{d(j)m) and has dimension I'km- The following well- 
known lemma, whose proof is recalled in [5] , gives a useful decomposition of E* in the hyperbolic 
region. 

Lemma 1.9. The stable subspace W{T,rj) admits the decomposition: 

(1.9) E^(r,7?) = 0^61 Ker L(d0™) , 

and each vector space in the decomposition (|1.9p admits a basis of real vectors. 



The next Lemma, also proved in [5], gives a useful decomposition of C^ and introduces projec- 
tors needed later for formulating and solving the profile equations. 

Lemma 1.10. The space C admits the decomposition: 

(1.10) C^ = 0^=1 Ker L(d</)„) 

and each vector space in (jl.lOp admits a basis of real vectors. If we let Pi, . . . ,Pm denote the 
projectors associated with the decomposition (ll.lOp . then there holds Im A^ (0) L(d0m) = Ker Pm 
for all m = 1, . . . ,M. 

For each ?n, G { 1 , . . . , M} we let 

denote a basis of ker L(d(/)m) consisting of real vectors. In section [3l we construct an approximate 
solution u'^ of (II. 2p of the form 

(1-11) ul{x) = ^^am,k{x,-^\ rm,k, 

m£lk=l ^ ^ 

where the am,k{x, 6m) are C^ functions decaying to zero as \9m\ — ^ oo, which describe the propaga- 
tion of pulses with group velocity v^, (see Proposition l3.6p . Observe that if one plugs the expression 
(jl.lip of u^ into P{eus,dx)u£, the terms of order 1/e vanish, leaving an 0(1) error, regardless of 
how the (Tm k are chosen. The interior profile equations satisfied by these functions are solvability 
conditions that permit this 0(1) error to be (at least partially) removed by a corrector that is 
sublinear (in fact bounded) as \6m\ -^ oo. Additional conditions on the profiles come, of course, 
from the boundary conditions. 

For use in the remainder of the introduction and later, we collect some notation here. 

Notations 1.11. (a) Let := ij^^ x M^, Qt ■= ^ n {-oo <t<T}, bQ := W^ x M^ b^r := 

bn n {-co <t<T}, and set ujt := 1+^^ n {-oo <t<T}. 

(b) For s >0 let H'^ = H^{bO,), the standard Sobolev space with norm {V{x' ,9o))s. 

(c) L^H' = L^{xd,H%bn)) with \U{x,do)\L2Hs = \U\o,s- 

(d) CH' = C{xd,H'{bn)) with \U{x,eo)\cH^ = sup^^>o|f/(.,Xd,.)|H= = \U\oo,s (note that 

(e) C°'*^(M^+^ X M) = {V{x\ Xd, 9o) G C (M+, C^'^ {R'^ x M, M^))} where C^' denotes the space 
of M times differentiable functions with derivatives up to the order M bounded. 



(f) Similarly, H^ = F^(60t) with norm {V)s,t and L'^H^ = L'^{xd,H^), CH^ = C{xd,H^) 
have norms |C/|o,s,T, |^|oo,s,T respectively. 

(g)When the domains of Xd and (x' , 9q) are clear, we sometimes use the self-explanatory notation 
C{xd, H%x', Oo)) or L^ixd, H^{x', 9o)). 

(h) For r > 0, [r] is the smallest integer > r. 

(i) Mo:=3d + 5, 

The main result of section [2] is the following theorem, which gives the existence of exact solutions 
to the singular system ()1.3p . or equivalently ()1.5p . and the original system ()1.2p on a time interval 
independent of the wavelength e: 

Theorem 1.12. Under Assumptions \l.i\ \1.^ \1.6[. \1.7\ consider the quasilinear boundary problem 
([12]), where G(x',^o) e H'+^{bn), s > [Mq + ^], satisfies 

Supp G C {t > 0} . 

There exist eq > 0, Tq > independent of e ^ (0, eo]? and a unique Ue{x,9o) € CHt n Li^H^ 
satisfying the singular problem (jl.Sp . so that 

x' ■ P 



is the unique C^ solution of (ll.2p on utq ■ 

Remark 1.13. The regularity requirement s > [Mq + — ^] in the above theorem is needed in order 
to apply the singular pseudodifferential calculus introduced in |7]. 

We can now state the main result of this paper. This theorem is actually a corollary of the 
result for singular systems given in Theorem 14. 161 

Theorem 1.14. Under the same assumptions as in Theorem \ 1.1^ there exists Tq > and func- 
tions crm,k{x,Om) G C^{il,TQ) Satisfying the leading order profile equations (j5.3p and defining an 
approximate solution u" as in (jl.lip such that 

lim n^ - Mg = m L°°{u;to), 
e— >0 

where Ue £ C^{ujto) is the unique exact solution of ()1.2p . In fact we obtain the rate of convergence 

^d 



u1\l°°{u:t ) < ^£2^1+5, where Mi :-- 



2^' 



Theorem 11.141 can be recast in a form where the pulses originate in initial data at t = and 
reflect off the boundary {x^ = 0}. This requires a discussion similar to that given in section 3.2 of 
[6] to justify the reduction of the initial boundary value problem with data prescribed at t = to 
a forward boundary problem (with data identically zero in t < 0), so we omit that discussion here. 

1.3 Profile equations 

In d space variables {x",Xd), consider the quasilinear problem equivalent to (jl.2p 

d-l 
ddUe + '^ Aj{eUe)djUe = F{eUe)Ue in Xrf > 

B{eue)ue = G{x' ,0q)\^^^^ on Xd = 
Ue = in i < 0, 

9 



where G(x', ^o) decays to like {Oq)~^ (for some A; > 2 to be specified later) as 16*01 — )• oo. For ease 
of exposition we will begin by considering the 3x3 case, which contains all the main difficulties. In 
section [SJ we describe the changes needed to treat the general case. We define the boundary phase 
(pQ := j3 ■ x' , the real eigenvalues ojm of A{I3), and the phases (pm '■= ^Po + ^mXd as in (jl.Sp . where 
/3 & Ti. We assume that the eigenvalues w^ are pairwise distinct. For the sake of clarity, we also 
assume that ui and ws are incoming (or causal) and UJ2 is outgoing. (The same kind of arguments 
would apply if two of the phases were outgoing and only one was incoming.) The corresponding 
right and left eigenvectors of the real matrix —iA{f3) are denoted rj and Ij, j = 1,2,3. 

Below we frequently suppress e-dependence in the notation. For functions U{x^9o,^d) and 
V(x,6'o,^d), define 

d-l d-l 

l{de„d^,) :=%+ j;/3,i,(0)a,„ =% + i(/3)ae„ and Z(a) := 9, + ^^1,(0) 9, , 
j=0 j=Q 

d-l 

M{U, de,V) := Y, (^J (d^j(O) • ^) deoV • 
i=o 

Formally looking for a corrected approximate solution of the form 

<(x) = [U'^ix,eo,^d) + i:^^ix,Oo,^d)]\g^^^^^^^^, 
we obtain interior profile equations 

£°: C{de„d^^)U'+L{d)U' + M{U'>,de,U^) = F{0)U\ 
and the boundary equation 

(1.13) 6° : B{0) ^/°U,=o,c,=o = G{x', 9o) . 

Consider the first equation in (11.12p . A function l/(^{x,6o,^d), taking values in M.^ and assumed 
to be C^ for the moment, can always be written 

U^ = 0-1(2;, 6*0, Cd)n + 0-2(2;, 6*0, Cd)?'2 + 5-3(2;, 6*0, Cd) Ca- 
using the matrix [ri r2 r^] to diagonalize A{/3), we find that the scalar ai must satisfy 

(% - iVi dg,) ai=0,i = l, 2, 3, in {(x, ^o, ^d) : ^o € M, ^d > 0}. 
This implies that the iTj's have the form 

^i{x,9o,^d) = cri{x,0o + uJiCd) for some ai{x,0i). 
Using ()1.13p . we find 

(1.14) 5(0) I ^ aiix', 0,60) rA = G{x' ,60) - B{0){a2{x' ,0,eo)r2), 

\i=l,3 J 
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Remark 1.15. 1. We expect the ai{x,6i) to decay polynomially to as \6i\ — )• oo. To prove this 
we must formulate and solve profile equations for the o"i's. For this we use an approach inspired by 
the formal constructions in |10j and ^17j . 

2. Instead of ai{x,9i) we shall sometimes write ai{x,9) with the understanding that is a 
placeholder for 6*0 + Wj S,d when it appears as an argument of Uj. 

To get transport equations for the iTj's, we consider (|1.12|) (e^): 

(1.15) C{de„d^^)U' = - (^L{d)U'> + M{U^,de,U'')) + F{0)U'> := H^,9o,^d)- 
The corrector U^ can be written as 

U^ =h{x,eo,Cd)ri + t2{x,eo,^d)r2 + h{x,eo,Cd)r3- 
Diagonalizing again we find that the tj's must satisfy 

(1.16) {d^^ - ujideo)ti{x,0o,U) = h ■ J^ ■■= J^i{x,Oo,^d), i = 1,2,3. 
The general solution to (|1.16|) is 

rS.d 

(1.17) ti{x,eo,S.d) =T*{x,eo+ujiS,d)+ J='i{x,eo+uJi{Cd- s),s)ds, 

Jo 

where t* is arbitrary. This can be rewritten 
tiix,9o,^d) 

(118) =Tiix,Go + ^i^d)+ Ti{x,9o + uji{^d- s),s)ds+ Ti{x,9o + ujj{Cd- s),s)ds 

Kd 

= Ti{x, 00 + uJiU) + / Ji(x, 6*0 + uJiiU - s), s) ds , 

J oo 



provided the integrals in (jl.lSp exist. 

We will need the following modification of a classical lemma due to Lax [13]. We refer to [6l 
Lemma 2.11] for the proof. 

Proposition 1.16. Let W{x,6Q,^d) = '^i=iWi{x,9o,^d)fi be any C^ function. Then 

3 3 

j=l 4=1 k^i 

where X^. is the characteristic vector fielc^ 

d-l 
3=0 

and V^ for k ^ i is the tangential vector field 

d-i 



V,::=Y,{kM0)rk)d,,. 



1=0 



^This vector field is a scalar multiple of dt + VXtiiv^'^i) ' ^x" that describes propagation at the group velocity 
Vi; see (fLS)) . 
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We see from (J1.15p that Ti{x,9o,S,d) has the form 
(1.19) Ti{x, 6*0, Cd) = -X4,Mi - ^ 4 dk degdh - ^ d^^ di de^Om + ^A^k-^ VJ. dk , 

k l^m k k^i 

'lOl 



(1.20) 



where we recall crp{x, 9o, £,d) = '^^(x, ^o + ^p(^d)- The coefficients in (I1.19P are defined b> 

d-l d-l 

i=0 j=0 

Thus, we compute 

Ti{x,eQ + uji{id- s),s) 

= -{X^^(Ji + c\ Oi dgai - e\ Oi){x, Oq + uji^d) 

- X] 4 '^k{x, 00 + WiCd + s{^k - Wj)) deak{x, Oq + ujiCd + s(wfe - ^i)) 

kj^i 

- X^ dl^ ai{x, 00 + a;i^d) de(Tm{x, Oq + Wj^^ + s{ujm - Wj)) 

- X] dJ,mT/(x,6'o + WiCrf + s(^/ - Wi))9eo-m(x,6'o + a;i^d + s(a;m -tJi)) 
+ X^(4 - ^fc) '7fc(a;, 6*0 + uJiU + s(wfc - iOi)) . 

k^i 

We look for functions ai{x,6) that decay at least at the rate {0)~'^ . So we assume now and 
verify later that they have this property. Then the integral 

rid 

(1.21) / F^{x,eQ + uJi{id-s),s)ds 

Jo 

is sublinear in {60, Cd) (a condition that must be satisfied by U^ if eU^ is to make sense as a 
corrector) if and only if the sum of the first three terms on the right in (jl.20p is 0. In that case the 
integral (|1.2ip is actually bounded, since the remaining terms in (jl.20p have good decay in s. This 
sublinearity condition gives the profile equations for the ctj's: 

Xtp.ai + c- o-j de^ai - e- cjj = 0, i = 1, 2, 3 

(1.22) {a^{x',0,eo),^ = l,3)=B (G(x',0o),^2(x',0, 0o)) , 

(Tj = in t < 0. 

where ;B is a well-determined linear function of its arguments whose existence is given by Lemma 
11.91 and the uniform stability assumption (the matrix [B(0)ri 5(0) rs] in (|1.14p is invertible). As 
expected from the general rule of thumb, pulses of different families do not interact at the leading 
order, meaning that the evolution equations for the amplitudes <Jj's are decoupled. In Proposition 
13.61 we show that system (I1.22P is uniquely solvable on some time interval [0,Ti], that (T2 = and 
that (Tj, i = 1, 3, decay at the rate {9)~'' for some A; > 2 to be determined. 



''We refer to section [5] for the general case. 
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Remark 1.17. The equations ()1.22p and our assumption that the ai decay at least at the rate 
{0)^"^ imply that the integrals in (jl.lSp all do exist. This argument is made more precise below. 

Next we introduce an averaging operator E and a solution operator Rqo that will be useful 
in the error analysis of the next paragraph. Motivated by the form of Ti in (jl.lOp . we make the 
following definition. 

Definition 1.18 (Type J- functions). Suppose 

3 

(1.23) F{x, Oo^id) = ^ Fi{x, Oo^id) n, 

i=l 

where each Fi has the form 

3 3 

(1.24) Fi{x,9o,^d) = J2fl{x,eo+uJk^d)+ Yl 9imix,Oo+uji^d)hl^{x,eo+uJm^d), 

k=l l<m=l 

where the functions fl,{x,9), gl^{x,9), h\^{x,9) are real-valued, C^ , and decay along with their 
first order partials at the rate 0{{9)~'^) uniformly with respect to x. We then say that F is of type 
J^. For such F define 

EF{x, 9o, Cd) ■■= Yl ir}^ f h- P{^^ ^0 + ^j {id -s),s) ds j r^ . 
Remark 1.19. 1.) For F as in pr23]) - (fL2il) . we have 

3 

(1.25) Y.F = YFi{x,9o + Uiid)ri, where Fi{x, 9) := f^x, 9) + gl{x, 9) hl{x, 9) . 

i=l 

2.) Observe that F as defined in (jl.lSp is of type F (hence the terminology), provided the fjj's 
have sufficiently regularity and decay in 9. In that case, we obtain 

3 

EF{x, 9o,id) = - ^ {X<p^cri + c- CTj de^ai - e\ CTj) r, , where Oi = ai{x, 9q + Wj ^d)- 

i=l 

Remark 1.20. The definition of E can be extended to more general functions. For example, if 

3 

i=l 

where the Fi{x,9) are arbitrary continuous functions, the limits that define EF exist and we have 
EiF = F. For another example, suppose F is of type F and satisfies EF = 0. Define 

(1.26) RooF(x,^o,ed):=^( / Vi(2;,0o + w*(ed-s),s)ds) r^ . 
Then the limits defining Roo-F" and ERoo-F" exist and we have ERoo-F" = 0. 
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Proposition 1.21. Suppose F is of type J- and satisfies Ei^ = 0. Then Roo-F is bounded and 
C{de,,d^,) RooF = Roo C{de,,d^,)F = F = {I - E)F. 

Proof. It just remains to show Yiao C{dQ^,d^^)F = F. This fohows by direct computation of the 
integrals defining Rqo >C(30(j,5^^)F and the fact that when EF = 0, we have Fi{x,9Q + u}i{^d — 

oo), oo) = 0. n 

The next Proposition summarizes what we have shown. 

Proposition 1.22. Let F(x,9Q,^d) be a function of type T. 

(a) Then the equation C{dgf^ , d^^)U = F has a solution bounded in (Oq, S^a) if o-nd only if^F = 0. 

(b) When EF = 0, every C^ solution bounded in {0o,Cd) has the form 

3 

U = 2^ Ti{x, 9q + ijji ^fj) ri + Roo-F with Ti{x, 6) ^ C and bounded. 

1=1 

Here WA = Y:l=i nix, Oq + Ui Cd) n and {I - E)U = RooF. 
(c)If U is of type T then 

EC{de„d^,)ll = Cide„d^,)EU = 0. 

Proof. Part (a) follows from the form of the general solution given in (J1.17p , and the fact that when 
a function F of type F satisfies EF = 0, the integrals 



/•oo 

/ Fi{x,eo+aJi{£,d- s),s)ds 
Jo 



are absolutely convergent. Part (b) follows from Remark 11.191 and ERqo-F = 0. Part (c) follows 
directly from Remark 11.191 D 

With the leading pulse profile 

U^{x,9Q,id) = cri{x,9o +ujie.d)ri + (J2{x,9q + ^2^0)^2 + cf?.{x,9q + uj^iid)^ , 

we can rewrite the profile system (jl.22p in a form that will be useful for the error analysis as follows: 

a)EU^ =U^ , 

h) E (L{d) U^ + M{U^, deJA^) - F(0) Z^°) = , 

c)5(0)Z^OU,=o,€.=o = G(x',^o), 
d) ^/° = in t < . 

These equations can also be obtained by applying the operator E to the equations ()1.12p . and have 
the common structure of weakly nonlinear geometric optics equations, see e.g. [201 chapters 7 and 

9]. 
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1.4 Error analysis 

We end this introduction with a sketch of the error analysis used to prove Theorem I4.16[ which 
yields Theorem 11.141 as an immediate consequence. The iteration schemes for the singular system 
(|A.6p and the profile equations (jl.27p are written side by side in (j4.14p . (j4.15p . For s largqlj and 
some To > 0, the proof of Theorem 11.121 produces a sequence of iterates U^{x,6o), bounded in the 
space E^ uniformly with respect to n and e, and such that 

lim [/" = Us in E^ uniformly with respect to e € (0, eo] ) 

where C/g is the solution of the singular system (jl.Sp . On the other hand the construction of profiles 
in Proposition 13.61 yields a sequence of profile iterates ly(^''^{x,9o,^d) bounded in £^ (see Definition 
14. ip and converging in £^~ to a solution U'^ of the leading profile equations (|1.27p . By Proposition 
14.31 this implies that the rapidly varying functions U£'^\x,9o) := U^''^{x,9o, ^) satisfy 

lim U^'" = U^ in E^^ uniformly with respect to e G (0, eo] • 

1 
Thus, in order to conclude \U^{x,6q) — Ue{x^6Q)\ps-a < C 8'^^'^+^ and thereby complete the proof 

of Theorem 11.141 it would suffice to show: 



Tn 



(1.28) There exists C such that for every n, |Z^°'" - U^\j^.-3 < Ce^Mi+s , 

The statement (|1.28p is proved by induction in section [H It is natural to try to apply the estimate 
of Proposition 14. 171 to the difference We'" — C/""*"^, but the problem is that for any given n, Ue'^ 
does not by itself provide a very good approximate solution to the boundary problem (j4.14p that 
defines U'^'^^ ■ Indeed, substitution ofUe'" into (|4.14p (a) yields an error, call it R^'^^ {x , 9q) , that 
is 0(1) in E^ . Since ^'^'"+1 satisfies (j4.15p . the main contribution to R'^'^^{x,9q) is given by 

7^"+l(x,eo,x°) ^^6^6 



(1.29) 71"+^ := (/- E) (^L(a^)^^'^'"+^ + M(Z^'^'",90,^^"'"+^) - F{0)U 

One would like to solve away the main error term in (jl.29p by using Proposition 11.211 and con- 
structing a corrector U^'^~^^{x,9o,^d) such that 

(1.30) £(ae„,%)^Yl■"+l = -7^"+^ 

and then use a corrected approximate solution of (I4.14|) (a) of the form 

While such a corrector is given explicitly by W^'"+^ = Roo(— '^""''^), it is not suitable for the error 
analysis because, although bounded, Ue'"" does not lie in any of the E^ spaces. 
To see the reason for this, note that U^'"^ has the form 

3 

(1.31) U^''\x, 9o, ^d) = Y. ^r(^, ^0 + uJi ^d) n. 

i=l 



^We take s > 1 + [Mo + ^] in Theorem ITT61 
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Observe that if / is a function that decays (say hke \s\ ^) as |s| — )• oo, the primitive J^ f{s)ds 
itself decays to zero as \0\ — t- oo if and only if / has moment zero ( /_^ f{s) ds = 0). Since neither 
^o,n+i j^Qj, ^Yie term M(Z^°'", 96)qW°'"+^) in ()1.29p has moment zerco, the definition of Rqo shows 
that this choice oilAe'"' generally cannot lie in any E^ space. We first try to remedy this problem 
using an idea inspired by an argument in [2]. We replace Z^*^'" (and similarly Z//'^'"^^) by a function 
Up'^ defined by functions cr" with vanishing first moments, where 

o-"p(x,m) := Xp("i)(T"(x,m), 0<p<l, 

and Xpi"^ is ^ low frequency cutoff function vanishing on a neighborhood of of size 0(|p|) and 
equal to one outside a slightly larger neighborhooq^J. We show the estimate 

These "moment zero approximations" (Definition I4.4p are the pulse analogues of the trigonometric 
polynomial approximations, which can be viewed as produced by high frequency cutoffs, used in 
the error analysis in the wavetrain case in [6l section 2.5]. With this change the contribution to 

7/l,n+l r 

Ue from 

-Roo(/-E) (L(a,.)<"+'-F(oX"^ 

lies in a suitable EJj^ space, but there is a problem due to "self-interaction terms" of the form 

coming from the M term in (ll.29p . which do not have moment zero. Thus, we replace these terms 
by i'^l^p^o'^^p )v ^^ ™ (|4.28p . The transversal interaction terms cr^pdga'^tj , i i^ j already yield 
contributions in an EJj^ space (Proposition 14.10]) . 

Using moment-zero approximations introduces errors that blow up as p — )• 0, of course, but 
taking p = e for an appropriate 6 > 0, one can hope to control these errors using the factor e 
in eU^'"'^^. Indeed, this works and by the process outlined above we obtain a corrector Up'^e 
which, though it does not solve away R^~^^, solves away "ah but 0{y/p+ m\+2 )" of ^e^^ ™ ^i})^ 
(see ()4.35p for more details) . Setting p = e^ and choosing the exponent b so that ^ = M\+i (so 
b = 2M +5 ) ' ^^ ^^^ ^^^^ *o apply the estimate of Proposition (14.17J) to conclude 

'd 



|^u,n+i _ u^+^\ ., < C7e^^7r+^ where Mi 

^0 



2+=" 



Remark 1.23 (Uniformly stable shocks). There is an analogue of Theorem 11.141 for uniformly 
stable shock waves perturbed by pulses. Uniform stability for the shock waves problem is an 
extension of Definition 11.51 and dates back to Majda [15]. The case of shocks perturbed by highly 
oscillatory wavetrains was studied in [6, section 3]. In that case there is a separate expansion for 
the oscillating shock front (a free boundary) 



(1.32) 



V'e(x') ~ o- xo + e ( X°(a^') + e X^ {x' 



1 (^> ^'i^') 



More precisely, we refer here to the moments of profiles like i7"{x, 9) or products of profiles that appear in these 
terms. 

^•^The cutoff renders harmless the small divisor that appears when one writes the Fourier transform of the ^-primitive 
of (j"p in terms of (T"p(a;, m). 
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in addition to an expansion for the solution on each side of the front. An important difference in 
the pulse case is that the term x^^i^') is absent in (jl.32p . and of course ^^(^'j^o) is now decaying 
instead of periodic in 9q. The expansions of the reflected waves on either side of the front are 
similar to (jl.lip . 

The singular shock problem has the same form as in the wavetrain case (see equations (3.39) of 
[6] and [22]), and the profile equations again take the form of equations (3.60) in [6], except that 
every occurrence of x^{^') is replaced by 0. The solution of the large system for the leading profiles 
is now considerably simpler than equations (3.68) of |6], since all the interaction integrals in that 
equation are now absent. This reflects the fact that pulses of different families do not interact at 
the leading order while wavetrains do. However, it is necessary to estimate interaction integrals in 
the error analysis. As in the pulse problem with fixed boundaries, one can in the shock problem 
obtain a rate of convergence of approximate solutions to exact solutions as e ^ 0. 

2 Exact solution of the singular problem 

The goal of this section is to prove Theorem 11.121 and solve the singular system (11. 5p . This is 
achieved, as in [24^ section 7], by solving the sequence of linear problems 

a) 5,.,C/ri + ^i^.(e[/;) fd,^ + ^2^) [/;+! = F(eC/;) U^ , 
(2.1) ^=° ^ ^ 

b)B{EunK-'X,=o = G{x',eo), 

c) U^+^ = in t < 0. 

As for the case of hyperbolic boundary value problems, that is without the singular parameter 1/e 
in the differential operator, see e.g. [3ll3], the solvability of each linear system (j2.ip relies on some 
a priori estimates. Our main focus here is the derivation of a priori estimates that are uniform with 
respect to the wavelength e. For the reasons detailed in the introduction of |24j . the appropriate 
functional setting in which one can derive uniform estimates is provided by the spaces E^ defined 
in (jl.6p . The main difficulty is to obtain L°° estimates uniform in e. These cannot be obtained 
simply from uniform L^(xrf, H^{x', 9q)) estimates since the estimate of dxJJe in terms of tangential 
derivatives provided by the system (12. ip blows up as e — t- 0. Much of the analysis in this section 
is similar to [24, sections 5 and 7], except that we use here the singular pseudodifferential calculus 
of Appendix El This introduces some modifications for the regularity assumptions in the results 
below. In the proofs of this section we shall often refer to [24] in order to keep the exposition as 
short as possible. 

2.1 Main estimate for the linearized singular problem 

We consider a linearized problem of the form 



rf-i 



/^i^So 



(2.2) ^=0 ^ ^ ^ 

b)BieV,)Ue\,,=o = 9e, 
c) C/j = in t < 0, 

where {Ve)eG{o,i] is a given family of functions, and {fe,9e) represent some source terms. Our first 
main result is the analogue of [241 Theorems 5.1 and 5.2] and proves unique solvability with a 
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uniform Lp' energy estimate for ()2.2p . The main point is to keep track of the regularity assumptions 
on the coefficients V^. 

Theorem 2.1. Let sq := [{d+ l)/2] + 1. There exists 5 > such that, for all K > 1, there exist 
some constants Jo{K) > 1 and Cq{K) > such that the following property holds: if the coefficients 
(^e)ee{o,i] iri (H^l) satisfy 

(2.3) leVeliacf^^-^ < 6 , |^e|cO,A/o(n) + lKlc(//''0(RdxR)) + k^a:^d^£li°°(f^) ^ -^ ' 

then for all T > 0, for all source terms fe G L'^{Qt), Qe G L?'{bQ.T) vanishing for t < 0, there exists 
a unique solution U^ € L?'{Q-t) to (j2.2p vanishing for t < 0, and this solution satisfies 

(2.4) |e-^*C/,|o,o,T + ^ (e-^*C/.U,=o)o,T < Co(K) (- |e-^*/.|oAT + ^ (e-^*5.)o,T') , 

V7 VT V7 / 

/or a// 7 > 70 (K). 

In Theorem 12. H the space C^'^^^iyt) denotes the space of functions v{x,6) such that for all 
^d ^ 0, v{-,Xd,-) is bounded on M*^ x M with all derivatives up to the order Mq bounded, and with 
all bounds that are uniform in Xd- The norm is defined by 

blcO."o(Q) := sup sup \\d^>^0v{-,Xd,-)\\L--{RdxR)- 
Xd>0 \a\<Mo 

For fixed x^, the {x' , ^)-regularity of a symbol enables us to use some of the symbolic calculus rules 
listed in Appendix lAl 

Proof. The first main step in the proof of Theorem l2.1l is to show a global in time a priori estimate. 
In other words, we consider a smooth function Us solution to (j2.2p . and wish to show the estimate 
()2.4p with T = +00. We begin with the following result. 

Theorem 2.2 (Kreiss, Metivier [121 I19j). There exists 6 > such that, if Bg denotes the closed 
ball of radius 6 in R , there exists an m x m matrix-valued function 

ReC^{Bs xR'^x (0,oo)), 

homogeneous of degree zero in {^', 7) and satisfying: 

(a)R{v,e,7) = R{v,e,7r; 

(b) there exist C > 0, c > such that for all (^,^',7).' 

(2.5) Riv,C',j) + CB*{v)B{v)>cT, 

(c) there exist finite sets of C°° matrices on Bs xR x (0,oo), denoted Ti, Hi, and Ei such that 

{^)Ke{R{v,i',^)A{v,i',^)) = Y,Tl{v.i\l) (^^^'(^^'^''^) SK^°?',7)) ^'*^'^'^''^)' 
(ii) Ti, Hi are homogeneous of degree zero in (^',7), Ei is homogeneous of degree one; 

(iii)Hi{v,e,7) = Hi{v,e,7), Ei{v,e,7) = E*{v,e,i); 

(iv) there exists c > such that 

Y,Ti{v,e,l)Ti*{v,C',l)>cI, Hi{v,C',7)>cI, Ei{v,C',7)>c{\e\+7)I. 
I 

The dimensions of Hi and Ei can vary with I . 
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The parameter 6 is fixed according to Theorem 12.21 We then define the following singular 
symmetrizer for the boundary value problem (j2.2p : 



7^,,-,:=Op^'^(i^(£K,e^7)), 

where singular pseudodifferential operators Op''''^(a) associated with a symbol a are defined in 
Appendix |Al We observe, as in \24\ remark 5.2] that our symmetrizer is not self-adjoint on L'^{Q). 
However, the remainder T^e,^ — 7^* is 0(1/7) as an operator on L^, uniformly in e. 

Under the regularity assumptions (j2.3p of Theorem 12. H the results given in Appendix |A] and 
the arguments in j241 pages 164-165] give the following properties for the symmetrizer TZ^^-y: 

{a)\ne,yW\o,o<C{K)\W\o,o, 
{b)\[d,„ne,,]W\o<C{K)\W\o, 

(c) Re ((7^,,^A,7 +Al^ns,^)W,W) > c{K)j\W\lo, 

(d) Re (7^,,^ W, W) + C{K) {B{eVe) W)l > c{K) {W)l , 

where Ae^^ denotes the operator 

-7io(eK) -Y.A,{eV,) (d^^ + ^^ 
j=o ^ ^ 

Let us focus for instance on property {<£) in ()2.6p . Since 7^* = 7^£^^-|- 0(1/7), Garding's inequality 
(Theorem lA.ip shows that it is sufficient to prove that the symbol -R(eV^, ^', 7) + O B*{eVfr) -B(eV^) 
is positive definite, and this property is given by (j2.5p . Other properties in (|2.6p are obtained by 
similar arguments (applying Propositions I A. 7l I A. 81 or I A.9P . see |241 pages 164-165] for more details. 
We perform the change of function C/g — )• e~^^Ue in (|2.2p . multiply (j2.2p a) by 'R-^'y ^-nd take 
the real part of the L?' scalar product with eT^^Us- The estimates (12. 6p yielco 

|e-^*f/,|o,o + ^ (e-^*f/eU,=o)o < C{K) (- |e-^*/.|o,o + ^ {^-^'oe)^ , 

V7 V7 V7 / 

for 7 sufficiently large, that is for all 7 > 70 (-fC). 

A similar uniform a priori estimate is valid for the dual problem (which satisfies the backward 
uniform Lopatinskii condition) . Then the arguments of [Sj H] , namely existence of a weak solution 
and "weak=strong" by tangential mollification, yields well-posedness of the boundary value problem 
(j2.2p . Localization in time is achieved as usual by showing a causality principle ("future does not 
affect the past"), which holds in our context since the constant C{K) in our energy estimate is 
independent of 7. D 

The uniform Lp' estimate (12. 4p enables us to show an estimate in the space E^ defined in (jl.6p . 
The result is similar to |241 Corollary 7.1] with a slight improvement with respect to the norm in 
which the source term /^ is estimated. 

Theorem 2.3. Let sq := [{d+ l)/2] -|- 1. There exists (5 > such that, for all K >1, there exist 
some constants 71 (i^) > 1 and Ci{K) > such that the following property holds: if the coefficients 
iye)eG{o,i] ^^ (!2.2p satisfy ()2.3p . then for all T > 0, for all source terms fs G !?• {H^ {hO-T)) , 



''The detailed computations can be found in [241 corollary 5.2], and are the singular analogue of [3l[4] 
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g,, E H^{bQ.T) vanishing for t < 0, there exists a unique solution Ue G H^{Qt) to (j2.2p vanishing 

for t < 0, and this solution satisfies 

(2.7) 

|e-^*C/.|oo,o,T + \e-^'Ue\o,i,T + -^ (e~^*f/.U,=o)i,T < C\iK) (- |e-T*/.|o,i,T + ^ {e-"ge)i,T) , 

for all 7 > 7i(i^). 

Proof. The regularity of the solution [/^ can be obtained by using the same arguments as in [U 
chapter 7] , that is by commuting the system (j2.2p with a mollified version of the Fourier multiplier 
of symbol (7^ + |i^'p + A;^)^'^. The argument shows that the dx' and dg derivatives of C/g are in 
L^, and (j2.2p then shows that the dx^ derivative of C4 also belongs to L^. We thus only show the 
estimate (j2.7p . 

1. L^ estimate of tangential derivatives. Commuting (j2.2p with a tangential derivative 
dtan £ {c^xQ) • • • ) f^i^d-i ) ^So}) we need to control the commutators 

Y,[M^Ve),dtan] (dx, + ^^) t^. = J] (dijl^K) " dtanVe) {s dx, + (3j dg,) U, . 

When multiplied by e"'''*, this source term is bounded in LP'{Q.t) by a constant times |e~'''*C/£|o,i,r 
and can therefore be absorbed from right to left by choosing 7 large. At this stage, we have 



(2.^ 



"^*C/.|o,i,T + -^ (e-^*C/.U,=o)i,T < Ci{K) {- |e-^*/.|o,i,T + ^ {e-^'ge)iA 



for all 7 large enough. 

2. L°°[L'^) estimate, part 1. We extend /^ and g^ beyond time T, which does not affect 
the solution Ue up to time T. Doing so, we just need to prove the L°°{L'^) estimate ()2.7p for 
T = +00. We consider a cut-off function x'^ i^i the extended singular calculus, that is a smooth 
function satisfying the conditions ()A.3P given in Appendix |Aj The L°°{L?') estimate is first proved 
on (1 — Xs(-O)) (6~'^*C4), where we let from now on Xs(-^) denote the Fourier multiplier whose 
symbol is 

Since | A: /3|/e is dominated by (7^ + IC'P)"'^" on the support of 1 — x*^, the same arguments as in [2 
page 173] yield 

1(1 - xl{D)) (e-^ *[/,)! 00,0 < C{K) (|e-^*/.|o,o + |e-^*f/,|o,i) 

(2.9) < C{K) {^ |e-^*/.|o,i + -^ (e-^*5.>i) ■ 

3. L°°{L'^) estimate, part 2. It remains to estimate \xt{D) (e~'^*C/e)|oo,0) which uses the fact 
that /3 is a hyperbolic frequency. More precisely, we can fix some parameters (5 > and ^2 > such 
that for all v in the ball of radius 6 and for all (z, rj) that are (52-close to /3, there holds 

Q{v,z,riy^ A{v,z,ri)Q{v,z,rj) = diag {Xi{v,z,r]), . . . , XN{v,z,r])) , 

for a suitable invertible matrix Q, and the Aj's satisfy 

{<^ — C7 if 7^1... p 
>C7, ifj=p+l,...,iV. 

20 



Moreover, Assumption 11.61 shows that the (square) matrix whose column vectors are 

B{v) Qi {v,z,r]),... , B{v) Qp{v, z, rj) , 

is invertible (here the Qj's denote the columns of Q). 

With the above notation, we can follow the proof of [243 Proposition 7.3], and write xti^) {^~^^Ue) 
under the form 

xf(Z?)(e-^*C/,) = roW, 

where, here and from now on, tq denotes a bounded operator on L^(il) whose operator norm is 
independent of e, 7, and where each component Wj of W satisfies a transport equation 

(2.10) d^.Wj - Xj{eVe,D,)Wj = ro{e-^'f,) + ro{e-^'Ue) . 

In (|2.10p . \j{eVe,Ds) denotes the singular pseudodifferential operator of symbol Xj{eVe,z,rj) (as 
described in Appendix [A|) . 

In the outgoing case {j = p + 1, . . . , N), we multiply (j2.10p by Wj, integrate from x^ to +00 
and apply Garding's inequality (Theorem I A. ip . obtaining 

r+00 
(W,(x,))2 + 7 / {Wj{y))ldy 

<C / {Wj{y))o{e-^'fe{y))ody + C / (W,(y))o (e-^*C/,(y))ody . 






The contribution of Wj on the right-hand side can be absorbed on the left by using Young's 
inequality, and Theorem 12.11 enables us to control the L^ norm of C/^. We thus get 

(2.11) sup |W,|L,o < C{K) [- |e-^*/.|g,o + A (e""*5.)g) • 

The estimates in the incoming case are similar, except that we integrate from to Xd- We thus 

sup |W,|L,o< sup (W,U,=o)g + C(i^) f^|e-^*/.|g,o + 4(e-^*5.)gV 



get 



Using the same arguments as in [24^ page 178], we can use the uniform Lopatinskii condition and 
write 

'Wi|,,=o\ /Wp+iU,=o\ 

I =r-o ; +ro(e~^*(7e) + -ro(e-^*C/eU,=o), 

from which we derive the estimate 

sup (W,U,=o)^ < C{K) sup |W,|L,o + C{K) (e-^*g,)g + ^ |e-^V.|g,o • 

i=l,...,p j=p+l,...,Af 7 

We combine the latter inequality with (j2.1ip . and recall Xs(-^) (6~^*^e) = ^0 W, so we get 

|x^(Z))(e-^U)|oo,o<C(i^) (^|e-^*/.|g,o + (e-^*5e>g) ■ 
Adding with (12. 9p and (12. 8p , we complete the proof of Theorem 12.31 D 
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2.2 Construction of the exact solution 

We use the iteration scheme (j2.1|) to solve the nonhnear system (|1.5|) . As usual, the convergence 
of the iteration scheme follows from the combination of two arguments: a uniform boundedness in 
a "high norm" (here in the space E^ given in (jl.6p ). and a contraction property in a "low norm" 
(here in E^). The estimate of a solution in E^ will be provided by Theorem 12.31 above, and we 
indicate below how we obtain the estimate of a solution in E'^, s G N. 

Proposition 2.4. Let sq := [{d + l)/2] + 1 and let k £ N. There exists 6 > such that, for all 
K > 1, there exist some constants jk{K) > 1 and Ck{K) > such that the following property holds: 
if the coefficients (K)eg(o,i] in ^3^ satisfy (H^D and belong to L"^ {H^+\bnT)) n L°°(F'=(60t)), 
then for all T > 0, for all source terms fg G L'^ {H^^^^ {bQ.T)) , 9e £ H''~^^{bQ,T) vanishing for t < 0, 
there exists a unique solution Ue G LF' {H^~^^ {bVLr)) n L'^ {H^ {bQ.T)) to (|2.2p vanishing for t < 0, 
and this solution satisfies 



(2.12) \e-^iUe\oo,k,T + \e- 



'^'U,\o,k+i,T + ^{e-^'U, 



1 



e\xa=o)k+l,T < Ck{K) { - |e "^ /e|o,fc+l,T 



1 

7i 



H ;= (e ^ 5e)fc+i,T + |^e|L°o(Tyi>°°{6QT)) 



-Ity 



elO,fc+l,T 



7 



+ 



-7*^ 



£\x,i=Q\k+l,T 



V7 



for all J > jkiK). 



Proof. The proof is like that of \2A\ Theorem 7.2]. One commutes (j2.2p with a tangential derivative 
d°' of order 1 < |a| < k, and applies (tangential) Gagliardo-Nirenberg inequalities. When the 
additional fast variable 9 lies in the torus M/Z, these inequalities are given in [2H Lemma 7.3], and 
we claim that the exact same inequalities are valid when the fast variable lies in M. 

When commuting ()2.2p with a tangential derivative 9" , one applies Theorem 12.31 and needs to 
control the commutator 



in the norm \e~'^'' ■ |o,i.r- The e factor in front of Vg cancels the singular 1/e factor and we obtain 
the estimate 



-7* 



MeVe) [d^^ + 



/?,( 



;d° 



U, 



0,1,T 



<C{K)\e-^'U,\o,k+i,T 

+ C{K) I C/e 1/^00(^1,0 



'mr)) 



"7*1/ I 



k+l,T ■ 



The |e ^^Ue\Q,k+i,T term on the right hand-side is absorbed by choosing 7 large enough, and we 
are left with ()2.12p . (Estimates on the boundary are similar.) D 

We can then deduce the main estimate in the space El^ defined in ()1.6p (the proof is the same 
as that of [241 Corollary 7.2] and is based on the choice T = I/7 in Proposition 12. 4p . 

Corollary 2.5. Let k > Mq + [^] and Ki,K2 > 1. Then there exist a constant C{Ki,K2) > 0, 
a parameter eQ{Ki,K2) G (0,1] and a time T{Ki,K2) > satisfying the following property: if 
T < T{Ki, K2), if the coefficients (V^)eg(o,ii in (12. 2p belong to Ef, and satisfy 



(2.13) 



|^£|oo,A:,r + |KUd=o|fc+l,T < Ki , ledx^VelL^f^n^,) < K2 , 
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and if e < eQ{Ki,K2), then for all source terms fe E L'^{H ~^^{bQ,T)), 9e ^ H ~^^{bQT) vanishing 
for t < 0, there exists a unique solution Ug € E^ to (|2.2p vanishing for t < 0, and this solution 
satisfies 

(2.14) \Ue\ooAT + \Ue\o,k+l,T + VT{Ue\c,,=o)k+l,T<Ck{Ki,K2)(T\fe\o,k+l,T + VT{ge)k+l,T) ■ 

The parameter eo in Corollary 12.51 is chosen so that (|2.13p implies \£Ve\L°°{nT) — ^ where 5 is 
as in Proposition 12.41 

We are now in a position to prove our main existence result for the singular system (jl.Sp . The 
norm in the space E^ is defined by 

\'"\eI^ •= l^|oo,fe,T + \v\o,k+l,T ■ 

Theorem 2.6. Let K > and let k > Mq + [-^]. Then there exists a constant K' > 0, a 
parameter £o{K) G (0, 1] and a time T{K) > satisfying the following property: the iteration (|2.ip 
with U^ = is well-defined for < T < T{K) and satisfies 

VnGN, ye<£o{K), \U^\e^ + \U^\^,=o\k+i,T < K , \ed^,U^\Lo.^Ur) < K' . 

Moreover, the sequence (C/") converges towards a function Us in EjT , uniformly with respect to 
£ G {0,£o{K)]. The limit Us belongs to Ef, and is a solution to (|1.5p . 

Proof. The constant K' is chosen such that, if \U^\^k < K, and if furthermore jV^I^;*: < K, then 
one has 



(2.15) 



d-i 

3=0 



<K' , 



independently of e G (0, 1]. Then the parameter eo is chosen as £o{K,K') given by Corollary 12.51 
The time 'T{K, K') > is chosen accordingly. Assuming that the induction assumption 

"^jl^n, ye<£o{K), \Ui\E,^ + {Ui\,,=o)k+i,T<K, \ed^,Ui\L^inr) < K' , 

holds (this is trivially true for n = 0) , we can apply the estimate (|2.14p of Corollary 12.51 to the 
system (j2.ip and obtain 

\U:^'\e^^ + Vr(C/r'U.=o>fc+i,T < Ck{K,K') (T\F{eU^)U^\o,k+i,T + VT{G)k+i.T 

< CkiK, K') (t C{K) + Vf {G)k+i,T 



Since {G)k+i,T tends to zero as T tends to zero, we can choose the time T small enough so that 
the induction assumption implies 

Ve < £o{K) , \U^+'\e, + {U^'+'\^,=o)k+i,T <K. 

Our choice of K' in ()2.15p implies that the induction assumption propagates from the rank n to 
the rank n + 1 because \£ dx^U^~^^\icc(^Q^-^ < K' . 

The convergence in EjT is obtained by showing a contraction estimate in Ej,, which is obtained 
by applying Theorem 12.31 We refer to |241 page 184] for the details. The limit Us of the iteration 
scheme ()2.ip is a solution to (jl.5p , which yields the regularity Us G E!^ (see [31 chapter 9] for similar 
arguments) . D 
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3 Construction of the leading pulse profiles 

Observe that we can solve the system (|1.22p by solving instead 

Xcp^ai + c\ Oi deOi -e\oi = ^, z = 1, 2, 3 
(3.1) {oiix', 0, 0), i = 1, 3) = e (G(x', 0), a2(x', 0, Q)) 

(T,- = in t < 0. 



where all occurrences of Qi or ^o firs now replaced by 0. To solve (|3.ip we use the iteration scheme 

(a) X^^<+i + c\ ara,<+i = e^ ar, i = 1, 2, 3 

(3.2) (6) (ar+i(x',O,0),. = 1,3) =i3(G(x',0),^ri(x',O,0)) , 
(c) CTf +^ = in t < 0. 

We will prove estimates for (|3.2|) in a class of Sobolev spaces weighted in Q. These weights are 
introduced in order to get an explicit decay rate in Q at infinity. 

Definition 3.1. For s G N and 7 > 1 define the spaces 

V := {a{x,e) E l2(E^+i x M) : {e,d^,dgfa E L^ for |/3| < s, and a = in t < o} , 
andr; :=e^*r% 
with respective norms 

(3.3) \a\s:= Y. 1^''^ ^f^ Sf a|i2(,,e) and |a|,,^ := le'^^aU. 

|/3| = |(/3i,/32,/33)|<s 

We will let H^ and H^ denote the usual Sobolev spaces with norms defined just as in (j3.3p but 
without the 6 weights. These spaces and those below have the obvious meanings when a is vector- 
valued. 

Remark 3.2. We have 

where "~" denotes an equivalence of norms with constants independent of 7 > 1. The second 
equivalence follows from 

The next proposition is helpful for estimating the commutators that arise when deriving F* 
estimates of solutions to the linearization of the profile system ()1.22p . Define 

A" := |a E L2(M^+i X M) : ^^^^ a E L^ for |/3i| < s, Sf^ a E L^ for |/32| < s, 

a^' aG L"^ for {fSsl < s, a = in t < o|, 

with 

\l3i\<s \l32\<s \l33\<s 

and let define Af, := e'''* A* with the norm \a\\s := |e~'''*a|A'' accordingly. 
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Proposition 3.3. The spaces F" and A* are equal, and the norms \a\s and \a\As are equivalent: 
there exists a constant C, such that 



(3.4) |a|As < \a\s < Cs \a 



s "A'* • 



Proof. Clearly, \a\\s < \a\s. The remaining inequality is proved by induction on s. The case s = 
is clear. The square |a|^ is a sum of terms 

(3.5) [{el^' 9f2 gft „) (^/3i Q^2 ^ft^) ^^ ^0 ^ 

where |/3| = |/3i| + |/32| + l/^sl ^ s. The terms with |/3| < s are dominated by C|a|^s by the induction 
assumption. 

Consider now a term like (j3.5p with |/3| = s > 0. Either 2|/3i| > s or 2|/32,/33| > s- Suppose 
2 1/32) /^sl ^ -s- Perform integrations by parts to obtain terms of the form 



(3.6) C 1 5°^g) a ■ {9, 8,, deT' adxde< Cs 19^,) a\l, +51(9, 8,, de^ oW^ , \ai 



s, i = 1,2. 



and other terms (where powers of 6 are differentiated) that can be estimated using the induction 
assumption. The second term on the right in (j3.6p can be absorbed by |a|^. Integrations by parts 
show that the first term on the right is dominated by the sum of a multiple of \a\\s and a term 
that can be absorbed by |a|^. 

The remaining case 2|/3i| > s is handled similarly. D 

Remark 3.4. 1) The spaces T^, F^, A*, and A^ all have obvious analogues when L^(M^+^ x M) is 
replaced by L?'{VLt) in the definitions, where we recall the notation 

Ot = {{x,9) G M^+^ X M : t < r}. 

The corresponding norms are denoted by adding the subscript T: \a\s^T-, \<As,'y,Ti I^|a= > etc. The 
equivalence (j3.4p continues to hold for the norms restricted to Ot; as can be seen by a standard 
argument using Seeley extensions [4!]. 

2) The analogous norms of functions of b{x' , 9) defined on M'^ x M or on bO^x are denoted with 
brackets: {h)s^T-, (^)s,7,T) etc. We denote the corresponding spaces by 6F|,, 6A^, etc. 

3) By Sobolev embedding it follows that if the functions aj appearing in (jl.lOp lie in Tip for 
s > —^ + 3, then J^ as in (jl.lSp is of type J^. Indeed, we then have, for example, 0| dg^aj S H^ 
for an index t > -^. 

4) More generally, if the functions /|, gj^, h\^ appearing in (jl.24p he in T^ for s > ^^ + 2, 
then F as in (jl.23p is of type F. 

We set |a|oo := \o\l°°(q.t) when the domain VLt makes no possible confusion, and we define 
Wrp'"^ := {a{x,9) : |a|i,oo := ^ \d"^0a\oo < oo} . 

|a|<l 

Estimates for the coupled systems. We can now state the main existence result for solutions 

(3.7) v°'"+H^,^) = K^^'^^^^^^) 

to the sequence of linear systems ()3.2p . 
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Proposition 3.5. Let T > 0, m> ^ + 1 and suppose that G G hVl^ and V^'" S T™ both vanish 
in t < 0. Then the system (|3.2p /las a unique solution V*^'"^^ G F^ vanishing in t < with 
(72 = 0- Moreover, there exist increasing functions, 7o(^) andC{K) of K := |V'''"|m,T such that 
for 7 > 7o(i^) we /ia?;e 

.0,n+l| , (V '""^ )m,7,r ^.^^ ( {G)m,'y,T \V '"|m,7,T 



(3.8) |V^''^+V,7,T + 7=^^"^ < C{K) '""'^"' + 

Proof. 1. L^ estimate. Anticipating the extra forcing terms that arise in the higher derivative 
estimates, we first prove an L^ a priori estimate in the case where a forcing term fi{x, 9) vanishing 
in t < is added to the right side of each interior equation in (I3.2J) . Setting F{x, 0) := (/i, /2, /s) 
we claim 

(3.9) I V°'"+i|o , T + ^^°'"^'^°'^'^ < C{K') f ^^i^ + i^^^^ + '^"'"'"•^•^ 'l , 

V^ V 7 \/7 7 / 

where K' := jV'^'^li^oo- The latter estimate is obtained by considering the weighted function 
g-7tyO,n+i^ and by performing straightforward energy estimates on ()3.2p (a). The traces of a"^^ , 
z = 1, 3, are directly estimated by using (j3.2p (b). 

2. Higher order estimates. We use again the system (]3.2p in its original form. Using the 
equivalence of norms established in Proposition 13.31 we first apply the L^ estimate to the problems 
satisfied by 6 ai, where k < m. The forcing term in this case is 

Clearly we may assume \9\ > 1. Applying <\2>.9\\ we can absorb the terms on the right involving 
cr"+"^ to obtain 

(3.10) 1^^ VO'"+%,,,T + ^^'^°'"^'^°'"'^ < C{K') f i%Ll + l^"'"'-'^-^ ') for 7 > lo{K'). 



Next we estimate 9^,0""'*"^ for \a\ < m. The forcing term in the problem satisfied by d^ify^^^ is 



:)Q 1 n+1 
i 



The commutator is a finite linear combination of terms of the form 
(3.11) (a^,Vr) {d'^Jdea^^^), |ai| + jaal = \a\, \ai\ > 1. 

We estimate these terms using the following two observations: 

A. Suppose mi + m2 > -5^, rui > 0. Then the product {a{x , 9) , b{x , 9)) — )• a • 6 is continuous 
from Hp X Hjp'' -^ H^. 

B. Suppose |qi| + \a2\ < \a\ < m. Then 

By A (with mi = m — \ai\, m2 = |ai| — 1) and B we have 

103. lip [0,7,7^ < C'krim.T kj"^ |m,7,T < Ci^ |(t"+ |m,7.r- 



26 



Applying (j3.9p and absorbing terms from the right, we obtain an estimate hke (j3.10p for d^,V^'^~^^ 
with C{K') replaced by C{K). 

The 6 derivatives dga-i, k < m, are estimated similarly. Derivatives involving dx^^ are estimated 
in the customary way using the tangential estimates and the fact that x^ = is noncharacteristic 
for X^^. 

3. Existence and uniqueness. This follows easily from the above estimates since the principal 
part of the system (13. 2p is given by three decoupled vector fields. One can therefore integrate along 
characteristics. Equations (|3.2p (a).(c) and the fact that X^^ is outgoing imply a2~^ =0. D 

Next we show convergence of the iterates V*^'" to a short time solution of the nonlinear profile 
equations (|3.ip . 

Proposition 3.6. Consider the profile equations ()3.ip . where G G hV^ , m > —^ + 1, and vanishes 
in t < 0. For some <Tq <T the system has a unique solution V^ G F^ with 02 = 0. 

Proof. 1. The iteration scheme ()3.2p defines a sequence (V'''") in F™. Fixing i^ > we claim that 
for T* > small enough, 

(3.12) |V°'"U,T* + (V°'")^,T* < K for ah n. 



Indeed, first observe that 



\u\ 



m,-y,T < Ci \u\m,T < C'2 e'*' |u|m,7,T, 



and fix 7 > 7o(-f^) such that ^y^ > 2C{K) Ci for 7o(i^) and C{K) as in Proposition 13. 5i Assuming 
()3.12p holds for n < no, we find that it holds for no + 1 after shrinking T* if necessary, using the 
estimate (j3.8p and the fact that G vanishes in t < 0. This new choice of T* works for all n. 

2. Convergence of the iterates in Tj, to some V*^ for a possibly smaller To > now follows 
from (j3.12p by applying (|3.8p when m = to the problem satisfied by (V'^'""'"^ — V'^'"). In view of 
()3.12p and a classical argument involving weak convergence and interpolation, we thereby obtain 
a solution V*^ G F^ with, in fact, a trace that lies in bV^ . This argument shows that the iterates 
V°'" converge to V° in F™-^ D 

4 Error analysis 

Next we carry out the error analysis sketched in section II. 4[ In section 14.11 we define and 
derive estimates for moment-zero approximations. In section IT2] we estimate interaction integrals 
involving both transversal and nontransversal interactions of pulses; these estimates are used later 
to estimate the first corrector Z^/ig. Finally, in section IT3l we complete the proof of Theorem 11.141 
by proving the stronger result. Theorem 14.161 

4.1 Moment-zero approximations to U^ 

When constructing a corrector to the leading term in the approximate solution we must take 
primitives in 9 of functions a{x,6) that decay to zero as |^| — )• oo. A difficulty is that such 
primitives do not necessarily decay to zero as \6\ — )• oo, and this prevents us from using those 
primitives directly in the error analysis. The failure of the primitive to decay manifests itself on 
the Fourier transform side as a small divisor problem. To get around this difficulty we work with 
the primitive of a moment-zero approximation to a, because such a primitive does have the desired 
decay. 

We will use the following spaces: 
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Definition 4.1. 1.) For s > 0, we recall the notation (jl.6p . that is E^ := {U G C{xd, H^{x' , 6*0)) Pi 
L'^{xd,H^ {x',9q))}. This space is equipped with the norm 

|f/(3;,^0)|i=;j:=|f/|oo,s,T + |f/|0,s+l,T. 

2.) LetS^ ■.= {U{x,9o,U) ■ Me^ - supg,>o l^(-, •,Cd)|i?j < 00}. 
Proposition 4.2. For s > (d + l)/2 ^/le spaces E^ and 8^ are Banach algebras. 

Proof. This is a consequence of the Sobolev embedding Theorem and the fact that L°°(6$7t) n 
H^{bQx) is a Banach algebra for s > 0. D 

The proofs of the following two propositions follow directly from the definitions. 

Proposition 4.3. (a) For s > 0, let a{x,6) G E^ and set d-{x,6o,S,d) '■= o-{x,9o + u}S,d), w G R. 
Then a £ £^ and 

Ic^l^j <C\a\^s+i. 
(b) For a G £^, set ae{x,9o) := a{x,0o, ^). Then 

Definition 4.4 (Moment-zero approximations). Let < p < 1, and let (j) G C°°(M) have supp 4> C 
{m : \m\ < 2} with (j) = I on {\m\ < 1}. Set(j)p{m) := <t>{^) andxp '■= ^ — <Pp- For a{x,6) G L'^{^t), 
define the moment zero approximation to a, ap{x,6) by 

(4.1) ap{x, m) := Xp{'m) (t(x, m), 
where the hat denotes the Fourier transform in 9. 

Proposition 4.5. For s > 1 suppose a{x, 9) G F^"^ , and define a {x,9o,^d) '■= <^{x,9o + (jj£,d)- Then 

a) \a - o-plfs < C \a\^s+2 ^, 

b) \dxd^ - dx^apl^s-i < C \a\^s+2 ^/p. 

Proof 1. Recall that a £T'^ ^ 9'^^ d^"" d^^a{x,9) e L'^{x,9) for \(3\ < s, which is also equivalent 
to d^ dt m^^ a{x,m) G L^{x,m) for \/3\ < s. It follows that 

(4.2) a G r^+2 => aix,m) G H^+^{x,m) C H\m,H'+\x)) C L°°{m,H'+\x)). 
2. We have 



\a — a. 

\a\+k<s+l 



\a\+k<s+l 

= Yj / / \d^ m^ a{x,m) (l)p{m)\'^ dx dm 

I I I I ^ 11 J\m\<2p J 
\a\+k<s+l ' '— '^ 

<C \a{x,m)\jjs+iMdm<C\a\l,s+2{2p), 

J\m\<2'p ^ 



where the last inequality uses (14.21) . The conclusion now follows from Propostion 14.31 

3. The proof of inequality b) in Proposition 14.51 is essentially the same. D 
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Proposition 4.6. Let a{x,9) G H^, s >0, and let CTp he a moment-zero approximation to a. We 
have 



a) Wp\h!^ < C\a\H^ , 

C 



(6) IfaGT^, then \ap\ri, < — \(j\r?^ 



P 
Proof. Part b) follows from ()4.ip . Indeed, for |/3| < s, 

C 



d^ 9f 2 ^/33 



&pix,m)\L'2^ < ^klr^, 



since |9m Xpl ^ C/p^^. Taking /3i = we similarly obtain part a). D 

Next we consider primitives of moment-zero approximations. 

Proposition 4.7. Let a{x,9) G T?p, s > 2+3- Let a*{x,6) be the unique primitive of ap in 9 that 
decays to zero as \9\ — )• cxd. Then a* G T^ with moment zero, and 

(Tri rs 

(6) \al\rs <C7^^. 

Proo/. 1. Since ap{x,9) G Tf,, s > f + 3, we have |crp(x,6')| < C {9)~'^ for all (x,0). The unique 
^-primitive of dp decaying to zero as \9\ — )• oo is thus 



cr*(a;,0) = — / ap{x,s)ds= / (Tp(a;,s)ds. 



Moreover, we have 



(4.4) 5e(T* = cr„ =^ im (T* = 3%^ = Xp iS", so a* = -r^ a. 

f f f iiji 

Since \m\ > p on the support of Xp, this gives 

\a{x,m)\ 



(4.5) k;(x,m)|<C' 



P 



and ()4.3p (a) follows directly from this. From (j4.4p we also obtain a*{x,0) = 0. 
2. The proof of ()4.3p (6) is almost the same, except that now one uses 



QS (Xp 



C 
< 



ps+l- 



Proposition 4.8. Let a{x,9) and t{x,9) belong to H^, s > ^^. Then 

(4.6) \crT- {aT)p\H^ <C\a\H^\T\H^y/P■ 



^ 
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Proof. With * denoting convolution in m we have 

|cJT-(fTr)p|^j ~ Y^ |9°m'=(<T*f)(x,m)(l-XpM)li2(x,m) 

|a|+fc<s+l 
m\<2p 



< C / \{cr*T){x,m)\Hs(^^)dm 

i \a{x,m - mi)\Hs(^^)\f{x,mi)\Hs(x)dmij dm 



< C , 

\m\<2p 



D 



Proposition 4.9. Let a{x,9) and t{x,6) belong to T^, s > 2 +3 and let (frr)* denote the unique 
primitive of {aT)p that decays to zero as \9\ — )• 00. Then 



{aT);\H^<C- 



IJ, I IJJJ, 



p 

Proof. Since T^ is a Banach algebra, Proposition 14.71 implies (cr)* G T^ with moment zero and 

|(ar)p|Hj<C . 

Since H^ is a Banach algebra, the result now follows from Proposition I4.6l fa). D 

4.2 Estimates of interaction integrals 

Pulses do not interact to produce resonances that affect the leading order profiles as in the 
wavetrain case. However, interaction integrals must be estimated carefully in order to do the error 
analysis. 

The following propositions will be used in the error analysis for estimating terms related to 
U^ as in (jl.30p . where the Ti appearing there are given by (jl.20p : in particular, we must estimate 
primitives of products of pulses. In some of the estimates below we must introduce moment-zero 
approximations to avoid errors that are too large to be useful in the error analysis. We begin with 
an estimate of "transversal interactions". 

Proposition 4.10. Let t be the smallest integer greater than | + 3 and let s > 0. Let ai{x,6), 
(72 (x, 9) belong to F^ n H^ and define 

(4.7) u{x,eo,^d)-= ai{x,9o+uj^d + as)a2{x,9o + uj£,d + s)ds, 

J 00 

where uj, a are real and a ^ {0, 1}. With Us{x,9o) := u{x,9o, ^) we have 

\ue\E^ < C{\ai\jjs+i |cj2|r^ + W2\h^+i kilr^)- 
uniformly for e G (0, 1] . 
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Proof. 1. For fixed Xd and e we first estimate 
(4.8) 



^d/s 



(Ti (T2 ds 



E 



H^{x',eo) \a\<s 



Xd/e 



d°, (cTi (72 ) ds 



+E 

L'^ix'fio) k<s 



f-Xd/e 
/ 9,^JaiC72)ds 



L2(x',eo) 



Here and below o"i, o"2 and their derivatives are evaluated at the points indicated in (|4.7|) with 
£,d = —, unless explicitly stated otherwise. 

2. To estimate As{xd) we consider for |ai| + \a2\ = \a\: 



^d/e 



a:,Vi«V2ds 



< 



L2(a;',0o) 



^d/e 



\d7^id:?a,\^,^^,^ds 



L^eo) 



< 



/oo 
(kl|L°°(x') W2\h={x') + Wi\h={x') W2\l-^(x')) ds 
-oo 



L2(eo) 



< Ai^sixd) +A2,s{xd), 



where we have used a Moser estimate in the x' variable. Setting z = 9o + to ^ + a s, we obtain 



Ai,s{xd) = C 



|cri(x,2;)|2,oo(x') 



oo 
oo 



e a a 



dz 



H''(x') 



L^ido) 



< C / \ai{x,z)\Lcx,(^x,^\a2\L2(e,H^(x'))dz 



C\a2\ 



LH 



H^{x')) / kl(a;,2)|L°°{a'') (^) TT 

J— oo \^/ 



dz 



(4.9) 



< C\a2\L2{g^H^{x'))Wl{x,z){z) \l°^{x,z) < C'|0'2|_ff|,(a;',6») kl|r« ) 



where the last inequality uses Remark 13.41 The estimate of ^2,^(3;^) is similar. 
3. Recalling the definition of the E^ norm and using 



|o"2|c(xd,ff|,(z',0)) < C\a2\H^xa,H;^(x',e)) < C\a2\jj^+i 



{x,d) 



we obtain from (14.9 



(4.10) 



\Ai^s{xd) A2,s{xd)\c{xa) + \Al,s+l{xd) A2,s+l{xd)\L^xa) ^ C'(kllH^+i kalr^ + k2|j:/^+i |cri|r^). 

4. To estimate Bs{xd) in (|4.8|) we consider for A:i + A;2 = k: 

/ f^0o '^1 ^00 ^2 ds = ± (dg^jO-i) (T2ds + (boundary terms). 

J 00 J 00 

Each boundary term has the form 5™Vi 9™^c72, ttt-i + m2 < k, where s is evaluated at Xd/e- We 
estimate such terms using Moser estimates as follows: 

('^•11) \^TQ^'^i'^9o^'^2\L'2{x',eo) < C'(kilL°°(x,0) k2lHs-i(x',eo) + kilH=-i(x',eo) k2lL°°(x,eo))- 
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For the integral term setting z = 9q + oj-^ + s, we have 



(4.12) 



/ id^^ai)a2ds <C \d^^ai (x,{eo + uj—){l - a) + az) a2{x,z)\L2^^,)di 



< c 






LHeo) 



£ 



LHdo) 



< C|(Ti|i2(a,/^j^s(5i)) \a2{x,z){z) \l^(^x,z) < C" kl U|,(x',6») k2|r^- 

From (jilT]) and (jiJ2]) . we obtain parahel to ()i30]) : 

l-Bs(a:d)|c(xd) + \Bs+i{xd)\L2(xa) ^ C'(kil/fj+i kalr^ + \^2\h^+^ ki|r^), 

completing the proof. D 

The previous estimate of transversal interactions did not require the use of moment-zero ap- 
proximations. However, nontransversal interactions of pulses can produce errors that are too big 
to be helpful in the error analysis. Thus, we are forced to use a moment-zero approximation in the 
next proposition. 

Proposition 4.11. Let a{x,9) and t{x,6) belong to F^, s > | +3. For a,ui gM., a ^ set 



Then 



fix,Oo,(,d) ■■= / {crT)p{x,6o+uj(,d + as)ds. 






f{x,eo,^) 



Proof. The integral equals a ^ {a t)*{x,6q + ^^i^ + a)) so the estimate follows first by applying 



Proposition 14.91 and then by applying Proposition 14.31 

Corollary 4.12. Let a{x,6), t{x,9), and uj,a be as in Proposition \4- 11\ and set 



D 



9{x,0o,U) ■= / {(TpTp)p{x,eo+ujU + as)ds. 



Then 



g{x,9o,—) 



j-is— 1 



V 



Proof. First apply Proposition 14.111 and then Proposition I4.6r a) . 

Proposition 4.13. For s > f + 3 let a{x,e) G H!^, T{x,e) G F^^^ With Lo,a £R, a y^ set 

rid 
h{x,6o,^d) ■=cr{x,eo + u}^d) / d0f^T{x,eo+uj(,d + as)ds. 



D 



Then 



h{x,eo,—] 



< C |(t| h" \t\ h" 
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Proof. The integral is equal to a ^ t(x, Oq + £,d{^ + o.)) so the estimate follows from the fact that 
E^ is a Banach algebra together with Proposition 14.31 D 

In the next Proposition we must use a moment-zero approximation since t{x,9) may not have 
moment zero. 

Proposition 4.14. For s > f + 3 let a{x, 9) G H^, r(x, 9) G Pf,. With uj,a gR, a ^ set 



Then 



,U) ■=deoCr{x,9o+u}£,d) / Tp{x,9o +ujU + as)ds. 

J oo 



j{x,9o,—) 



<C — ^ 



E^-^ p 



Proof. The integral is equal to a Tpi^, 9q + ^d(a; + a)). The estimate follows by the argument of 
Proposition 14.13] except that now we also need Proposition 14.7( a) and Proposition 14.6( a) . D 

The proof of the next Proposition is evident from the proof of Proposition 14.141 

Proposition 4.15. For s > | + 3 and uj,a gR, a ^ 0, let a G Fy and set 

r<d 



Then 

Xd 



rKd 
k{x,9o,^d)= ap{x,9o + u}^d + as)ds. 

J oo 

k{x, 



^5 ^U; 



\C\ H" 

4"' P 



4.3 Proof of Theorem [TTl 

Now we are ready to prove Theorem 11.141 which shows that the approximate solution u'^{x) 
converges in L°° to the exact solution n^ of Theorem 11.121 as e — )■ 0. In this section we prove the 
following more precise Theorem, which implies Theorem 1 1 . 141 as an immediate corollary. As before 
we focus on the 3x3 strictly hyperbolic case to ease the exposition. The mostly minor changes 
needed to treat N x N systems satisfying Assumptions [TTTl 11.21 and 11.61 are described in section [5j 

Theorem 4.16. For Mq = 3d + 5 and s > 1 + [Mq + ^], let G{x', 9o) G bVp-^ and suppose G = 
in t < 0. Let Us{x,9o) G E^ he the exact solution to the singular system ()1.5p for < e < Cq 
given by Theorem \1.1S\ let V^ = ((Ti,(T2,(T3) G T^ be the profile given by Proposition \3^ and let 
U^ G £^^ be defined by 

3 

^° (2;, 6*0 , Cd) := X] '^J (^' ^0 + Wj id) Tj . 
i=i 

Here < Tq < T is the minimum of the existence times for the quasilinear problems (|1.5p and 
(fOT]) . Define 

U^{x,9o):=U''{x,9o,^). 
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The family U^ is uniformly bounded in E^ for < e < Sq; moreover, there exists < Ti < Tq and 
C > such that 



(4.13) \U,-U^\^s-3 < Ce^Afi+s^ 

where Mi is the smallest integer > -^ -\-3. 

The proof of Theorem 14.161 will use the strategy of simultaneous Picard iteration first used by 
Joly, Metivier, and Rauch in [llj to justify leading term expansions for initial value problems on 
domains without boundary. Consider the iteration schemes for the quasilinear problems (|1.5p and 



(4.14) 



and 



a 

j=0 

c) U^+^ = in t < 0, 



/0,n^ 



(4.15) 



a)EZ^°''^+i=Z^°'"+i 

6) E (l(9)Z^°''^+i + M(Z^°'", 99,^^°-'^+^)) = E (F(0)ZY° 

c)S(O)ZYO'^+^U,=o,c,=o = G(x',0o) 
d) U^'^'^^ = in t < 0, 

where U^'"'{x, 6q, ^d) '■= Yli=i '^7i^-: ^0 + ^j ^d) fj for a" as constructed in Proposition 13. 5i Setting 

ZYO'"(x,0o):=Z^°'"(^,^o,y), 

we observe that to prove the theorem it suffices to prove boundedness of the family lA'^ in E^ along 
with the following three statements: 



(a) lim [/" = Ue in E'|, uniformly with respect to e S (0, eo] 



(4.16) 

K,^ rr" — TT ;,, i^ 

(b) lim W°'" = U^ in i?^"-*^ uniformly with respect to e G (0, eo] 

n— >-oo 

(c) There exist positive constants Ti < Tq and Ci, independent of n, such that for every n 

The first statement, together with uniform boundedness of the families U^, Us in E^ , is proved 
in Theorem 12.61 by showing convergence of the scheme ()4.14p using the following linear estimate 
(which is a consequence of Proposition 12. 4p . 

Proposition 4.17. Let s > [Mq + — 5— ] and consider the problem (|4.14p . where G G H'^^^ vanishes 
in t <0, and where the right side of (|4.14p i^a) is replaced by J- & E^. Suppose C/" G E^ and that 
for some K > 0, ei > 0, we have 

\K\e^ + I e £>,.,[/; I LOO <K foree (0,ei]. 
Then there exist constants Tq{K) and eo{K) < £1 such that for < e < eo o.iT'd T <Tq we have 
|C/;+ib. +Vt{U^+X^=o)s+i,t < C{K) {t\T\e^^ + Vt{G)s+i,t) ■ 
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Proof of Theorem \4-16[ 1. The boundedness of Z^/^ in £^|, and (j4.16p (b) follow directly from Propo- 
sition liT3| together with the fact that V° G H^^'^ and V"'" -^ V° in H^^. (In fact, the proof of 
Proposition 13.61 shows that the V*^'" are bounded in r^"*" and V*^'" — )• V'^ in F^ .) 

2. The approximate solution Ug'^ is by itself too crude; in order to prove (j4.13p using Propo- 
sition [5T7] we must construct a corrector eZ^/ig that lies in some EJp space. To achieve this we 
first approximate U^'"^ and ^o,n+i j-^y nioment-zero approximations Up'^ and Up'^ . For now we 
fix < p < 1 and define for each n 

3 

3=1 

where a'^ is the moment-zero approximation to a"^ defined by (|4.ip . Thus we have Up'"'{x, Oq,^^) = 
'EUp'^{x,9o,S,d) and by Proposition 14.51 

(4.17) |WO'"-Z^°'"|^j-i<CvP, and \d,,l('''''+' - d,,U^'"+\s^^^2 < C ^ , 



for C independent of n^i. 

3. For now we express the induction assumption as: there exists < a < 1 (to be determined) 
and positive constants Ci, Ti < Tq such that 

(4.18) \U^^-U^'''\j^s-3<Cie''. 

The boundedness of the family U^ in E^ together with ()4.18p imply 

\F{eU^) U^ - F(0)Z^°'"|^s-3 < Ce". 

To 

In view of (j4.17p and Proposition 14.31 this implies 

(4.19) \F{eU:)U:-F{ml'^\E^---^<C{VP + n- 

4. Define 

Gp := Z(a,)Z^O'"+i + M(^/0'",5,„Z^°'"+i). 

We claim that 

(4.20) \^Gp - E(F(0)Z^°'") 1^.-2 < C^. 

To 

Indeed, from (j4.17p and the explicit formula (jl.25p for the action of E on functions of type J-, we 
have 

|E (F{<d)U^''' - F(0)Z^°'") Ls-i < CVP. 

To 

But E(F(0)Z^°'") is given by the left side of (lil^ fb). so ^^?M follows by observing that (fiTT]) 
and Proposition 14.21 implv 



(4.21) 



E (L{d^) (^0'"+i -ZY0'"+1)) < C^ 

Tq 



'Constants C, Ci... appearing in this proof are all independent of n and e. 
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Here we have used the fact that the arguments of E in (j4.2ip are functions of type J-", so the formula 
(jl.25p can be appHed. 

5. Next define the operator 

1 



Lo := L(9,) + ^L(d</)o)9eo + M(ZY°;,", ^^o), 



which is an approximation to the operator appearing on the left side of (j4.14p (a) that will allow us 
to use Proposition 11.211 to construct a useful corrector lA}^. Indeed, we claim 



(4.22) 



.^U:+^ - F{eU^)U^\j,s-z < C(vP + 

^0 



This follows from ()4.14p (a) and the estimates 



\MeU^)d,^U^+^ - i,(0)9.^.[/;+i|^.-i < Ce, 



(4.23) 



-^MeU^)(3jdg,Ur' 



- A,iO) (3, de.U^-^' + dA,{0) ■ U^ (3, de,U, 



n.+ l 

e 



<Ce, 



T7>S — 1 

^0 



dA,(0) • {U: - W^'^) /3, de,Ur' ^._3 < C \U^ - <'rUr3 <C{^ + e«). 



0,n+l 



6. Construction of the corrector. First observe that since £{80^^ , d^^)Up' = 0, we have 



and thus 

(4.24) 

We have 

(4.25) 



^0 Uple — Sp,e > 



Lo<r+i - F{0)U^:: = Gp,, - F{0)U'^:: 



{EiGp - Fimin), + {{I - ^Gp - Fimin), ■ 



E(gp-F(0)^0'"))J^^-. <Cv^, 



by (j4.20p . the formula (jl.25p for E, and Proposition 14. 3i The second term on the right in (j4.24p 
is not small, so we construct U}^ to solve (most of) it away. By Proposition 11.211 the function 



U, 



1 ._ 



0,nN 



Roo (/ - E){Qp - F{<d)Up'^) satisfies 



(4.26) 



C{de,,d^,)Ul = -{I- E)(gp - F{Q)U^^ 



0,n\ 



However, this choice of Up is too large to be useful in the error analysis. 

7. To remedy this problem we replace (/ — E)^p by a modification [(/ — E)Gp\mod defined as 
follows. First, using (J1.19p and Remark 11.191 we have 

(4.27) (/ _ E)g, = ^ [ - j; 14Vgi + Y. A^lpde^'^Z' + E ^U^P^So^-y 

where a^p = agp{x,9Q + ujqid)- The problem is caused by the nontransversal interaction terms 
given by the middle sum over k ^ i, so we define 



(4.28) [(/ - E)gpUd = E - E ^^<^ + E 4(^^,p5eo<+^)p + E dl,m<rr,pde,aZ+p' 

i=l \ k^i k^i ly^m 



r-,;, 
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and we set 

(4.29) U^ := -Roo ([(/ - E)gpUod -{I- E)F(0)^/°'")) . 
Instead of (j4.26p we have 

(4.30) C{de,,d^,)Ul = -[{!- ^)Gp]mod + {I - E)F(0)^/0'". 
For later use we set 

D{X, eo,^d) ■■= {I - ^)Qp - [{I - ^)Qp\mod 

and estimate 

(4.31) |D(x,0o,y)ls-3<Cv^. 
Indeed, using Propositions 14.81 and I4.6r a) we have 

I (^.V^o<^' - {<pde,all\) {x, 9o + a;,.^)|^^-3 

8. Estimate of \U^Jj^s-2. By g22|) , (g^S]) and the formula (fT^ for Roo, for i = 1,2,3 we 

must estimate \bi{x,6Q,^)\^a-2, where bi{x,9o,^d) 



^T 



(4.32) 

ri. 



5Z 4 / ' (<p5eo<p') (a:, 00 + ^iid + ^(wfc - Wi)) ds+ 

2_/ "^Im / '^*>(^' ^0 + ^iS.d)d0oCrl^pix, 9o + LOiid + s{ujm - Wi)) ds+ 

Joo 

zl '^Irn / cr"p(2;, ^0 + ^iCd + s(w/ - a;i))50(,(Tj^+^(x, 6*0 + uji^d + s{u)m - ^i)) ds+ 

XI ^k / crfc,p(3^' ^0 + ^iCd + s{uJk - Wi)) ds- 

X] / ^fc^fc,p(^' ^0 + ^iCd + ^(Wfc - Wj)) ds = X bi^rix, 6*0, Cd), 



l^i •'"^ 



rid 



k^i-"^ r=l 

where 6^^^, r = 1, . . . , 6 are defined by the respective lines of (|4.32p . Since V*^'" is bounded in H^ , 
using Corollary 14.121 we find 

\fjn\ Ifl ^"+^1 

Ki(x,go,^)U.-.<CE '"'''' [' ""'""^ <C/P. 

k^i ^ 

Similarly, from Propositions 14.131 and 14.141 we get respectively 

\bi,2ix,0o,—)\E^'^ <C, \bi^^{x,eo,—)\^s-2 <C/p. 
e T e T 
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Since V '" is actually bounded in Tp' , Proposition 14. lOl on transversal interactions implies 



e -"T 
where we have used Proposition 14.6( b) to estimate 



C 



Err — rt+1 



P'- 



(t:. 



n+l| 
m,p Ir, 



t+i < 



C 






T p'- 

By Proposition 14.151 we have \bi^^{x,9Q,^)\^s-2 < C, and the estimate of bi^ is the same, so 
adding up we obtain 

C 

2 < - 



(4.33) 



K 



p,elEtr^ - ^t+1 



p- 



To estimate {dxJJ}))e we differentiate (|4.32|) and estimate as above to find 



\{dxJjl)e\ES^A 



< 



C 



pt. 



(4.34) 

9. We claim 
(4.35) |Lo {U'^i^+'+eU'p^,) - F{0)U^;^\^s-j <C (^ 

Indeed, we have 



p' 



,t+2 



so by ^^2^ and ()i30]) we find 

(E(g, - F{0)U^'^))^ + D{x, 00, ^) + {meU},)e + M{U^;2, de,){eU^,,). 



The estimate ([1:22]) now follows from (lOSD . (lOTT) . (103]) . and I^M). 
Using (|iT9]) . (HJH, and ([05]) . we obtain 

(4.36) |Lo (C/r^ - (^°;r 1 + 5ZY,\J) |^._3 <C{^ + e^ + ^ 

10. Next we claim that the following estimates hold: 

(4.37) 

/3 5eo, 



)• 



(a) 



(^9., + A(eZ^O;;, a., + ^)) ([/;+! - (^^0;;+i + e^^ij) 



<C(vP + e" + -T-2) 



"To 



P' 



(6) |S(£Z^°:r) (C/r^ - (Z^°:r+ ^ + eU'^J) |^.-. < C (VP + e" + ^^). 

To P 



Indeed, (|4.37p (a) follows from (|4.36p by estimates similar to (|4.23p . while (j4.37p (b) is a simple 
consequence of (|4.14p (b) and (|4.15p (c). Applying Proposition 14.171 we find 



\Ur' - «r' + eUl,J^s^ <C^o{Vf^ + e^ + ^), 



p- 



and thus 



\U. 



n+l _ 7/0,n+l 



-■i<C^/To{^ + e'' + 



ii+2' 



Recall that t is fixed and equals Mi in the notation of Theorem 14.161 Setting p = e we compute 



,^^ ziT5' ^° ^^ lake a = 2 = 2^ ^^"^ complete the induction step by shrinking 

To to a small enough Ti if necessary. This completes the proof of Theorem 14.161 D 



Vp = W^ when b - 2j^, 
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5 Extension to the general N x N case 

Here we describe the relatively minor changes needed to treat N x N systems satisfying As- 
sumptions ll.H II. 2^ and II. 6[ We first describe the construction of profiles in the general N x N 
case. For each m £ {1, . . . , M}, let 



"m,k 



, A; — 1 , ... , i^fc^ , 



denote a basis of real vectors for the left eigenspace of the real matrix iA{f3) associated to the real 
eigenvalue —uJm and chosen to satisfy 

{1 , if m = ?n,' and k = k' , 
U, otherwise. 

For V € C^ we set 

Pm,k V := {im,k ' v) Tm^k (no complex Conjugation here). 

Functions of type J- (see Definition ll.lSp have the form 

(5.1) F{x, ^0, Cd) = X] X] ^^M^^ Go,Cd) rm,k 

m=l k=l 

where each scalar function Fm^k is decomposed as 

(5.2) Fm,k = ^ fZ' (^' ^0 + ^m' U) 

m' 

+ Yl 5'm',fc',m",fc" (^^^0+ ^m' Cd) K^i ^k' ,m" ,k" {^^^0+ ^m" id) ■ 



m'k'Tn",k" 



In (j5.2p . m' S {1, . . . , M}, k' G {1, . . . ,Vk /}) and similarly for (m", k")\ moreover, the functions 
fm' k' ^^^- have the same properties as the corresponding functions in Definition ll.181 The averaging 
operator E is given by 

EF := V] ( lim — / F,n,k{x,0o + u>rniU - s),s) ds ] rm,k , 
^ V^^- T Jo J 

and for F as in (|5.ip . it follows that EF = J2m k Fm,k fm,k where 

Fjn,k ■■= fm^{x,9Q+U}mid)+ ^^ V'fc'.m.fc" (^' ^0 + ^m Cd) /i^;t,m,fe" (^' ^0 + ^^ ^d). 

k',k" 

On functions of type J- such that EF = 0, the action of the operator Rqo is given by 

R-ooF ■■=Y[ I Fm,k{x, 00 + UJmiid - s),s)ds\ rm,k ■ 
m,k ^-^"^ ^ 

The general form of the profile equations (|1.27p still applies. With 

W{x,eo,id) = y^^Wm,k{x,9o,id)rm,k, 

m,k 
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the decomposition of Proposition 11.161 now has the form 

L{d)W = y^^{X^^Wm,k) rm,k + X! yZ ^Z''!k''^rn',k' \ rm,k, 

m,k m,k \m,' ^m,k' 

where V^'^, is the tangential vector field 

rf-i 
^m'\k' '■= XI (^'"'*^ ^^(^^ r^',k') dxy 

3=0 

In place of (|1.3ip and (|3.7p we now have 

M fkm 

m=l k=l 

^ ' J m=l,...,M;k=l,...,Uk^ 

The argument that led to the profile system ()3.2p now giveqlj 
(5.3) 

j=Ok,k'=l k=l 

(c) 0-^+1 = in t < for ah m, k, 
where the coefficients b^ ^ . are defined by 

(5.4) b^ij := £rn,i ■ /3j (dij(O) rm,k) rm,k'- 

Remark 5.1. There is a potentially serious obstacle to proving estimates for the system ()5.3p . If 
one takes the L^ pairing of (|5.3p (a) with a^i {x, 9), it is not clear how to use integration by parts 
in 9 to move the 0— derivative in the sum on the left onto the n-th iterate. This problem does not 
arise in the estimate for p.2p . The next Proposition, which is p, Proposition 2.18], removes this 
difficulty by showing that there is a symmetry in the coefficients that appears after regrouping. 

Definition 5.2. For u near let —ujm{u), m= 1, . . . ,M, he the eigenvalues of 



d-l 

' 'u) 






and Pm{u) the corresponding projectors. 

The functions 00m{u) and Pm{u) are C°° for u near since j3 then belongs to the hyperbolic 
region of A{u.,^'). 



^^The nonlinear equations for the functions am,i are, of course, obtained from (|5.3I) by removing the superscripts 
n and n + 1. 
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Proposition 5.3. Let w £ M be expanded as w = ^^ j^ Wm,kfm,k = Xlm ^"i o-nd define 

d-l l^km 

(5.5) 5z>W:=EEC^™.'=' 

j=0 fc=l 

where the h^ ^ ■ are defined in (|5.4p . Then there holds 

t.R\ vim, N /-da;m(0)-u;™ ifk' = l, 

(5-6) 5;fc,(w;) = <^ , 

lU otherwise. 

Proof. We differentiate tlie equation 

d-l \ 

Wm(u) / + E ^i ^i ("") -Pm(n) = , 
i=o / 

witti respect to u in the direction Wm, evaiuate at n = 0, and apply Pm '■= -Pm(O) on tlie left to 
obtain 

d-l 

(5.7) Pm E ^J (d^i(O) • '^ni) Pm = i-dcOmiO) ■ Wm) Pm ■ 

j=0 



The second equality in (j5.5p and (j5.7p imply (j5.6p . D 

Proposition 15.31 allows us to write 

d-l l^kn^ 

(5-8) J2 E C '^-.^ ^^^:S' = B^ii^'n 5^<y ' 

i=0 A;,fc'=l 

where W"'" := Ylmk'^mk^rn,k] hence we can shift the 0— derivative and integrate by parts as 
discussed in Remark |5. 11 Using (j5.8p . we deduce from (j5.3p that o"^"*], = when m S C Otherwise 
the proof of Proposition 13.51 goes through as before. The statement of Proposition 13.61 is thus 
unchanged, except that in the second sentence we have c^^k = when m € O now. 
The formulation of Theorem 14.161 is exactly as before except now 



IVl i^km 

U^{x,9o,U) = E yZ^rn,k{x,9o + u)mU)rm,k, V^{x,9) = (ajn,k{x,9)] 

^ — ■ ^ — ' \ / m=l, 



m=l k=l 

The error analysis in the proof of Theorem 14. 161 goes through with the obvious minor changes. For 

'■:,p'^9o"k,p 



example, the troublesome self-interaction terms c\ c^„ Seo'^fc p , /c / i, in (I4.27P are now replaced by 



terms of the form c^ ^ ^, o"^ ^ ^do'^^k' n' "^ t^ ^' where the index p as before denotes a moment-zero 
approximation. These terms are handled just as before by introducing [(/ — 'Ei)Q]mod-, see ()4.28p . 
in which they are replaced by cl^/^f., i^^m k p^So^^k' p^v '^^^ contribution of these terms to the 
corrector Up,, is estimated as before using Corollary 14.121 
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A Singular pseudodifFerential calculus for pulses 

Here we summarize the parts of the singular pulse calculus constructed in j7] that are needed 
in this article. First we define the singular Sobolev spaces used to describe mapping properties. 

The variable in M'^"''^ is denoted {x,9), x E M'^, G M, and the associated frequency is denoted 
{S,,k). In this context, the singular Sobolev spaces are defined as follows. We consider a vector 
/? G R"' \ {0}. Then for s G M and e G ]0, 1], the anisotropic Sobolev space ii'^''^(M'^+i) is defined by 

and / (1+^+— I \u{^,k)fd^dk<+oo}. 



Here u denotes the Fourier transform of u on M^+^. The space H'^'^{W^'^^) is equipped with the 
family of norms 



V7>1, y u € H^'^R'^^') , ll^lll^.,,^^ := ^^-^ / (72 + 



-1 



e 



kP ^\ 



\u{i,k)?didk. 



When m is an integer, the space //"''^(M +^) coincides with the space of functions u G L'^iMr^^] 
such that the derivatives, in the sense of distributions. 



(8..f*) 



"1 / fl \ "d 

dxi + -de] u, Qi H h Od < m , 



belong to L'^{M.'^^^). In the definition of the norm || • ||//m,e ,|,, one power of 7 counts as much as one 
derivative. 

A.l Symbols 

In this Appendix, O denotes an open set and nolonger denotes the set of outgoing phases. Our 
singular symbols are built from the following sets of classical symbols. 



Definition A.l. Let O C M be an open subset that contains the origin. For m ^ R, we let S™(C') 
denote the class of all functions a : O xM.'^ x [1, 00) — )- C*^, M > 1, such that a is C°° on O xM.'^ 
and for all compact sets K C O: 

sup sup sup(7' + |C|')-('"-'^l)/'|5,-5|a(t-,e,7)l <Ga,.,K- 

vGK fGiRd 7>1 



Let C^{W^^^), A; G N, denote the space of continuous and bounded functions on M'^+^, whose 
derivatives up to order k are continuous and bounded. Let us next define the singular symbols. 

Definition A. 2 (Singular symbols). Fix /3 G M'^ \ {0}, let m £R and let n G N. Then we let 5^ 
denote the set of families of functions (a£^-y)£gio,ii,'y>i that are constructed as follows: 

(A.l) V(x,0,e,fc) gM'^+^ xM'^+S ae,^ix,9,C,k)=a(eV{x,e),C + —,7 

where a G S"^{0), V belongs to the space C^(]R^^) and where furthermore V takes its values in a 
convex compact subset K of O that contains the origin (for instance K can be a closed ball centered 
round the origin). 
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All results below extend to the case where in place of a function V that is independent of e, 
the representation (jA.ip is considered with a function V^ that is indexed by e, provided that we 
assume that all functions e V^ take values in a fixed convex compact subset K of O that contains 
the origin, and {Vi;)s£{o,i] is a bounded family of C^(M'^"''^). Such singular symbols with a function 
Vg are exactly the kind of symbols that we manipulated in the construction of exact solutions to 
the singular system (jl.3p . 



A. 2 Definition of operators and action on Sobolev spaces 

To each symbol a = (ae,7)e6]o,i],7>i S 5^ given by the formula (jA.ip and with values in C^^^, we 
associate a singular pseudodifferential operator Op^'"'{a), with e G]0, 1] and 7 > 1, whose action 
on a function u G 5(IR^^; C ) is defined by 

(A.2) Op^'^(a)n(x,e):=— i-^ / e'^^'^-^'^'^ a (eV{x,e),^ + ^,^) u{(,k)d^dk . 

Let us briefly note that for the Fourier multiplier a{v,^,j) = i^i, the corresponding singular 
operator is dxi + (/3i/e) dg. We now describe the action of singular pseudodifferential operators on 
Sobolev spaces. 

Proposition A. 3. Let n > d + 1, and let a G S^ with m < 0. Then Op^''^{a) in l\A.2\\ defines a 
bounded operator on L'^(M.'^^^): there exists a constant C > 0, that only depends on a and V in the 
representation (jA.ip . such that for all e G ]0, 1] and for all j > 1, there holds 

Vug 5(M'^+1) , llOp^'^(a) u||o < -^ ||n||o . 

The constant C in Proposition IA.3I depends uniformly on the compact set in which V takes its 
values and on the norm of y in C^^ . For operators defined by symbols of order tti > we have: 

Proposition A. 4. Let n > d + 1, and let a G S™ with m > 0. Then Op^''^(a) in ()A.2p defines a 
bounded operator from H"^''^ [W^~^^) to L'^{M.'^^^): there exists a constant C > 0, that only depends 
on a and V in the representation (jA.ip . such that for all e G ]0, 1] and for all 7 > 1, there holds 

VnG5(M'^+^), ||Op^'^(a)u||o < C ||u||h™.S7 • 

The next proposition describes the smoothing effect of operators of order —1. 

Proposition A. 5. Let n > d + 2, and let a G S^^ . Then Op'^''''(a) in <\A.2\i defines a bounded 
operator from L^(IR^^) to H^''^{W^^^): there exists a constant C > 0, that only depends on a and 

V in the representation (jA.ip . such that for all e G ]0, 1] and for all j > 1, there holds 

VnGcS(M'^+^), ||Op"'^(a)'u||^i.,^^< C||n||o. 

Remark A. 6. In applications of the pulse calculus, we verify the hypothesis that for V as in (jA.ip . 

V G C^(M'^+i), by showing V G H'{M.'^+'^) for some s> ^ + n. 

A. 3 Adjoints and products 

For proofs of the following results we refer to [7]. The two first results deal with adjoints of singular 
pseudodifferential operators while the last two deal with products. 
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Proposition A. 7. Let a = a{eV,S, + ^,7) G 5°, n > 2 (d + 1), where V G i/''o(M'^+^) for some 
sq > -^ — h 1, and let a* denote the conjugate transpose of the symbol a. Then Op^''^(a) and 
Op^''^(a*) act boundedly on L^ and there exists a constant C > such that for all e g]0, 1] and for 
all 'J > 1, there holds 

Vue5(M'^+M, ||Op^'T(a)*w-Op^'T(a*)u|L < -||n||o. 

7 

If n> 3d + 3, then for another constant C , there holds 

Vug 5(M'^+1) , ||Op"'^(a)* u - Op^'^(a*) uW^,,.^^ < C \\u\\o , 

uniformly in e and 7. 

Proposition A. 8. Let a = a{£V,C + ^,-f) e S^, n > 3d + 4, where V G F^o(]R'^+i) for some 
So > -^ — h 1, and let a* denote the conjugate transpose of the symbol a. Then Op'''"' (a) and 
Op^''^(a*) map H^'^ into L^ and there exists a family of operators i?^''^ that satisfies 

• there exists a constant C > such that for all e G ]0, 1] and for all 7 > 1, there holds 

VmGcS(M"'+^), \\R^"Iu\\q<C\\u\\q, 

• the following duality property holds 

Vu, t; G 5(R^+^) , (Op"'^(a) u, v)l2 - {u, Op='^(a*) v)l2 = (i?"'^ u, v)l2 . 

In particular, the adjoint Op^''^(a)* for the L^ scalar product maps if^'^ into L^. 

Proposition A. 9. (a) Let a,b £ S^^, n > 2(d + l), and suppose b = cr{£V,S, + -^,j) where 
V G II'^"{M.'^^^) for some sq > -^ + 1. Then there exists a constant C > such that for all 

£ G ]0, 1] and for all 7 > 1, there holds 

Vug 5(M'^+^) , ||Op^'^(a) Op^'^(6) u - Op^'^(a6) u\\q < - \\u\\o . 

If n > 3d + 3, then for another constant C, there holds 

Vug 5(M"'+^) , ||Op^'^(a) Op^'^(6) u - Op^'^(a b) uW^i.e^^ < C \\u\\o , 

uniformly in e and 7. 

(b) Let a G 5"^,^ ^ S^ or a & 5"^, 6 G 5"^, n > Sd + 4, and in each case suppose b = 
cr(eF, ^ + ^,7) where V G H^°{W^~^^) for some sq > ^^ + 1. Then there exists a constant 
C > such that for a// e G ]0, 1] and for all 'j > 1, there holds 

Vug cS(M^+1) , ||Op^'^(a) Op^'^(&) u - Op^'^(a6) u||o < C \\u\\o ■ 

Our final result is Garding's inequality. 

Theorem A.l. Let a & S^ satisfy Kea{v,^,^) > Ck > for all v in a compact subset K of O. 
Let now a ^ Sq, n>2d + 2 be given by (|A.lj) . where V G -H'*o(M'^+^) for some sq > ^ + 1 and is 
valued in a convex compact subset K . Then for all 5 > Q, there exists 70 which depends uniformly 
on V , the constant Ck o-nd 6, such that for all 7 > 70 and all u G 5(M'^+"'^), there holds 

Re (Op^'^(o) u; u)l2 > {Cr - 5) \\u\\l . 
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A. 4 Extended calculus 

In our proof of L°^{X(i;L^{x',6q)) estimates for the linearized singular system (Theorem | 
we use a slight extension of the singular calculus. For given parameters < 6i < 62 < 1, we choose 
a cutoff x^(C', ^,7) such that 

< X' < 1 , 

(A.3) x^(e',^,7)^ion{(,^ + ie'iY/^<.,|M|}, 

suppx^c{(7^ + |e?)^/^<<^2|^|}, 

and define a corresponding Fourier multiplier Xd ^^ the extended calculus by the formula (|A.2p 
with X*^(C') jT) ill place of a{£V,X,^). Composition laws involving such operators are proved 
in [7], but here we need only the fact that part (a) of Proposition IA.9I holds when either a or 6 is 
replaced by an extended cutoff x*^- 
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